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Abstract 


This paper is a critical survey of the measurement of rainfall by remote 
sensors. One parameter radar rainfall measurement is limited because both 
reflectivity and rain rate are dependent on at least two par^eters of the 
drop size distribution (DSD), i.e., representative raindrop size and number 
concentration. This paper extends earlier work to develop a generalized 
rain parameter diagram which includes a third distribution parameter, the 
breadth of the DSD, to better specify rain rate and all possible remote 
variables. Simulations show the improvement in accuracy attainable through 
the use of combinations of two and three remote measurables. 

The spectrum of remote measurables is then reviewed. These include path 
integrated techniques of radiometry and of microwave and optical 
attenuation. One carefully designed short path microwave attenuation 
experiment is sufficiently persuasive to show that the disappointing 
results achieved in many others was due largely to a combination of rain 
sampling problems and vertical air motions between the path and the gages. 
Also, when paths are col inear, attenuation deduced from radar and 
radiometry are in good agreement with that measured directly. By and 
large, dual wavelength radar methods which were aimed at improved 
measurements in small range increments have proven disappointing. However, 
when the attenuation estimated in this way, or by radiometry, is used as a 
constraint on the retrieval of rain profiles from the radar, the results 
are more promising. Selected experiments involving combinations of two or 
more of the three measurables, radar reflectivity, attenuation, and/or 
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radiometry. show considerable promise when adequate account is taken of the 
sampling and air motion problems. The best results in gate-by-gat'-e 
measurements have been achieved with dual polarization or differential 
reflectivity However, even these have failed to meet their full 

potential because rainfall often does not behave according to apriori 
assumptions. An accompanying paper shows that the use of a third remote 
parameter in addition to Z and offers great promise. In many cases, 
accuracy can be great ''v enhanced and ambiguities resolved by personal 
examination of the raaar displays which depict the nature of the bright 
band, cell spacing, homogeneity, and drop size sorting effects. With a 
growing appreciation of the needs and the capabilities of the various 
techniques, the future of highly improved remote rainfall measurements 


seems bright. 


1. Introduction 


The multi-parameter measurement of precipitation is the result of the 
relatively recent recognition of the facts pointed out many years ago that 
both the rainfall rate, R, and the radar reflectivity factor, Z, were 
functions of at least two parameters of the drop size distribution (DSD), 
namely representative drop size and nimiber concentration. Indeed, both R 
and Z are also functions of the breadth of the DSD. As a result, it is not 
possible to utilize a single Z-R regression relationship for all types of 
rain; nor is it possible to account for space and time variations from the 
mean Z-R law which are due to fluctuations in the DSD induced by a variety 
of physical processes. This paper reviews the work done to measure 
rainfall more accurately without any apriori assumptions about the DSD 
through the use of additional remotely measurable parameters such as 
microwave and/or optical attenuation (Ulbrich and Atlas, 1978), 
differential reflectivity (Seliga and Bringi, 1976), and microwave 
radiometry . 

We concentrate almost exclusively on rainfall and exclude snow and 
hail except as the characteristics of the former may be used indirectly to 
infer features of the rain. We allude to hail mainly when it may be 
confused with rain. We also note that this is a selective, not a 
comprehensi ve review. Thus, while we shall refer to the original work for 
basic concepts, we shall omit many others and cite only those wtiich appear 
to have been pivotal in the development of a branch of thought. 


Background 

It is over 40 years since Ryde 19^6) first showed the 

relationships between the reflectivity and attenuation of microwaves and 
rainfall rate and thus set the foundations for radar measurements of 
precipitation, and indeed, for all of radar meteorology. It was Wexler and 
Swingle (1947), Marshall et al . (1947), Atlas (1947, 1948), and then 
Marshall and Palmer (1948) who turned Hyde's work, which was concerned with 
precipitation as unwanted clutter, toward the meteorological applications 
of radar. Following Marshall et al. (19^7), Austin and Williams (1951) 
were the first to conduct thorough quantitative experiments on the relation 
between echo intensity and rain rate. Austin and Richardson (1952) then 
noted that the relations depended on the rainfall type. This, and the 
subsequent work by Blanchard (1953) in which he showed major differences 
between orographic and non-orographic rains in Hawaii were among the first 
signs that there was not a universal relationship between the reflectivity 
factor, Z, and rain rate, R. 

In spite of the persuasive evidence which had already been gathered by 
the mid-50*s (Twomey, 1953) that one could not hope to determine R from a 
measurement of Z alone and a simple Z-R relation, much of the community 
persisted in this vain hope. Of course, it was also hoped that one could 
find Z-R relations appropriate to either the regional climate, to the 
observable weather parameters, or to rainfall type. Attempts to devise 
such classification schemes have been reported by Cataneo (1969) and by 
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Cataneo and Stout (1968). The result is that by 1970 Battan (1973) 

reported 69 Z-R relations from many places around the world. It is safe to 

say that that number has probably doubled in the last decade. Further 
attempts to resolve these problems have involved averaging in space and 
time the rainfall rates determined from Z-R relations to deduce storm 
totals. Wilson and Brandes (1979) review this work and conclude that 
considerable averaging is required to obtain radar measured rainfall which 
is within a factor of two of the true rainfall 75 % of the time. 

In recognition of the wide variation in Z-R relations, Atlas and 
Chmela (1957) proposed a physical basis for those relations and developed a 
rain parameter diagram (see also Atlas, 1964) involving the four basic 
parameters, Z, R, W (rain water content), and (the median volume drop 

diameter). Plots of the various observed Z-R laws on the diagram showed 

dramatically that the variations were due largely to differences in the 
association between and R and that the latter did not usually correspond 
to that found by Marshall and Palmer (1948) and Atlas (1953) for stratiform 
rainfall. In short, when rain is comprised of many small drops of low 
fallspeed, the liquid water content is bound to be higher and the 
reflectivity lower than with an equal rain rate of large fast falling 
drops. The result would be an erroneous value of R as deduced from a Z-R 
relation. Moreover, these variations can occur from day to day, place to 
place, and moment to moment, especially as a result of size sorting. 

Among other things, Atlas and Chrnoia U957) showed that rain rate, R, 
could be expressed as a function of Z, and G, where the latter quantity 



a dimensionless moment of the mass distribution which increases with the 
''vreadth of the distribution (Atlas, 196^). In other words, an accurate 
sptcif ieacion of R depends not only upon the measurement of Z but also upon 
the measurement of both and G, or proxies thereof. 

Thesn facts should be recognized as well from the nature of any 
drop size distribution (DSD) sucn as that of Marshall-Palmer, viz., 

N - exp (- AD) ( 1 ) 

where and A are parameters of the distribution with A expressible in 

termz of D through AD = 3.67 (Atlas, 1953 ). Marshall and Palmer ( 19 ^ 8 ) 
o ^ o 

take to have a specific value independent of rainfall rate, viz., 

4-3-1 -0 21 

= 8 x 10 m cm and A to vary as A = 41 R * or, alternatively, 

0 21 

D^ = 0.09 R * . Embodied in their result is the explicit assumption that 

one of the distribution parameters is constant so that the distribution is 
a function of only one variable parameter, viz., D^. If this assumption 
were strictly true then all integral quantities defined in terms of N(D) 
would be dependent on only one parameter and would therefore all would be 
uniquely expressible in terms of any other integral parameter. That such 
is not the case is clear from the work of Waldvogel (1974) and Donnadieu 
(19H0) who show that and D^ can undergo sudden, independent, very large 
changes in rainstorms. Consequently, these integral quantities are 
functions of two independent parameters and possibly more if it is 


necessary to specify distribution parameters such as the breadth and the 
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upper and lower limits and respectively. 

Having recognized that the remote measurement of rainfall rate 

req lired at least two independent parameters, Ulbrich and Atlas ( 1978) 

generalized the diagram of Atlas and Chmela (1957) to include additional 

parameters as shown in Fig. la. This diagram contains the isopleths of Z, 

R, W, and D as in Atlas and Chmela (1957) but also carries isopleths of N 
O o 

and optical extinction . On separate overlays isopleths are shown for 
total number concentration and microwave attenuation A ()) at radar 
wavelengths ^ = 0.86, 1.25, 1.778, and 3.22 cm. The temperature 
dependence of A(^) is also included. A typical attenuation overlay is 
shown in Fig. 1b for A = 3.22 cm. Since Seliga and Bringi (1976) showed 
that Zpp, the differential reflectivity factor at orthogonal polarizations, 
is a unique function of for an M-P DSD, the isopleths of in this rain 
parameter diagram may be replaced by isopleths ^DR using the inset in 
Fig. la. 

llie nomogram of Ulbrich and Atlas is a powerful tool. Most 
importantly, by entering the diagram with any two of the 8 parameters (or 9 
if is included), the other 6 (or 7) are uniquely determined without 
apriori assumptions other than the basic one that the distribution is 
exponential. (See below for a generalization to other distribution forms,) 

The diagram also demonstrates that an assumption of an empirical 
relationship between any two rainfall parameters necessarily implies 
relationships between all other pairs of parameters on the diagram. This 
conclusion is substantiated by Ulbrich and Atlas (1975) who show that all 
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of the empirical relations which can be established betv#een pairs of 
rainfall parameters involve scatter in the experimental data about these 
relations which is directly related and of about the same relative 
magnitude. For example, the percentage deviation of an experimental data 
point from an empirically derived Z-R relation will be of about the same 
magnitude and have the same sense as the deviations of this data point from 
the corresponding empirically derived - R relation. 

An illustration of the use of the diagram is shown in Mg. Ic where 
several empirical Z-R relations have been plotted on a simplified version 
of the rain parameter diagram. Also shown as a shaded region is that area 
within which lie all of the 69 Z-R relations listed by Battan (1973). As 
indicated by Ulbrich and Atlas (1978) this shaded region spans a range of 
values of almost 3 orders of magnitude and, for given Z, a range of R 
values of more than one order of magnitude. In addition, these Z-R 
relations generally cross both the and isopleths indicating that both 
Z and R are dependent on two distribution parameters each of which varies 
independently between successive observations. Attempts to improve the 
accuracy of prediction of Z-R relations by adjusting the coefficients and 
exponents will therefore be largely unsuccessful. Deviations found from 
use of such equations are the result of the natural variability in both N 

o 

and and cannot be simultaneously accounted for in a Z-R relation, but 
can be allowed for in dual-measurement methods. 

From the above, it is apparent that we need to measure at least one 
independent parameter in addition to Z to obtain more accurate measures of 
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rainfall and to do so without any aprlori Z-R relation. In a series of 
papers, Ulbrich and Atlas (1975, 1977) and Ulbrich (1981), have simulated 
the increased accuracy attainable through the use of combinations of 
microwave attenuation. A, and Z, optical extinction, and Z, and in an 
accompanying article in this volume (Ulbrich and Atlas, 1982a), the 
combination of differential reflectivity Z^^^ and Z. All of these will be 
reviewed in Section 4 where it will also be shown that for some 
dual-measurement methods the determination of two parameters of the DSD 
results in a residual variance and/or bias in the specification of rain 
rate. Thus, in order to attain the “ultimate" in accuracy, one should 
measure a third parameter which is a function of the shape or breadth of 
the DSD. The effects of allowing for deviations in DSD shape from 
exponential on the structure of the rain parameter diagram are covered in 
Section 3. In Section 4, it is shown that by introducing a third remote 
measurable in experimental simulations, it is possible to eliminate the 
aforementioned systematic bias due to DSD shape variations. 

The simulations of multi-p'^rameter measurements by Ulbrich and Atlas 
are all based upon the computation of the remote measurables from several 
hundred samples of DSD's. Thus they do not include the measurement errors 
and practical limitations of actual remote sensing systems. Consequently, 
the "ultimate" accuracy demonstrated by these simulations assumes that the 
measurement technique involves combinations of perfectly measured remote 
measurables and not those which would be derived from a very dense raingage 
network which presumably woul i provide the ultimate measure of ground 
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truth. ^ We cover these problems In suounary fashion in Section 5. 

Also Included in Section 5 is a sunmary of the various remote methods 
of measuring rainfall which have been proposed or actually utilised in the 
field. Included in this section is a brief history of single measurement 
techniques involving radar reflectivity, microwave attenuation, or similar 
measurables. The results obtained thus far from dual-measurement methods 
are also reviewed, which includes dual-wavelength radar methods, 
dual-polarization methods, and combinations of radar and radiometer 

measurements. Special attention is paid to the apparent success of methods 
»A)ich involve path-averages. Excluded from this summary are direct methods 
such as raingages or even vertically pointing Doppler rads.' located within 
the rain region as a means of calibrating the remote sensors or bounding 
the errors in correction algorithms. The former have been discussed by 
Hitschfeld and Borden (195^)) and by Wilson and Brandes ^1979). However, we 
shall discuss several approaches by vdtich we may be able to infer features 
of the rainfall which are not directly measurable in Section 6. 

3. Generalization of Rain Parameter Diagram 

In the next section it is shown that the assumption of a two parameter 
exponential size distribution in a dual-measurement method of determining 

^The sampling errors associated with individual gages and the gate spacing 
relative to storm size raise serious questions as to what really 
constitutes "ultimate" ground truth. 
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rainfall rate Is, in most cases, adequate to achieve at least a five-fold 
increase in accuracy over that which results from the use of an empirical 
relation. However, for some combinations of remote measurables further 
Improvement is possible if allowance is made for deviations of the 
experimental size spectrum from the assumed exponential form. The 
implication is that to achieve the maximum accuracy of measurement requires 
knowledge of the shape of the size distribution. Ulbrich and Atlas (1982a) 
describe methods by which the shape of the distribution could be determined 
either apriori or by direct measurement. They also present evidence for 
varying size distribution shape in nature through analysis of empirical Z-R 
relations of other workers. Further support for such shape effects is 
contained in the work of Blanchard (1953), Best (1950), and Mueller (1965). 
The latter work Is especially significant since it involves rain drop size 
spectra collected by a photographic method that does not require knowledge 
of the drop fallspeed to find the size distribution and, in addition, has a 
much larger sampling volume than other methods. This method is used to 
find average drop size distributions from several locations around the 
world which in every instance are better described by a mathematical form 
which is closer to a gamma or log normal distribution than to exponential. 

To show the effects of varying distribution shape on the structure of 
the rain parmneter diagram it may be assumed that the distribution is of 
the form 


N(D) = N d" exp (-AD) (0 < D < <«>) 
o — 


( 2 ) 
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where the exponent m may have any positive or negative value. Such a form 

means that the distribution is now described by three parameters (N^, ^ and 

m) and that one of them must be known apriori in a dual-measurement method 

or that a third measurable would have to be introduced if all three 

parameters were to be determined simultaneously. It should be noted that 

for this size distribution the parameter A is related to D^, the median 

volume diameter, by the excellent approximation ® = (AD ) = 3.67 + m. 

m o m 

The import •'nt moment G introduced by Bartnoff and Atlas (1951) is given by 

which varies monotonically from 5.2 to 1.8 as m varies from -2 to 2, 
respectively. As indicated by Ulbrich and Atlas (1982a), negative values 
of m correspond to broad, concave upward distributions, such as those found 
in orographically induced rainfall in Hawaii, whereas positive m 
corresponds to a narrow, concave downward distribution such as those found 
by Mueller (1965) for a variety of rainfall types and by Jones (1956) for 
thunderstorm rain. 

In terms of the size distribution N(D), the quantities plotted on the 
rain parameter diagram are defined by equations specified in detail in 
Ulbrich and Atlas (1978). The definitions of W, Z, E , and all yield 
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results which can be expressed simply as a product of the three factors N , 

o 

a power of D^, and a function of m when Eq. (2) is used for the size 
distribution. The integrals defining R and A would be done numerically 
unless it is sufficiently accurate to adopt power law approximations of the 
form v(D) = Y D°*^^(Yr 17.67m s"’cm“°*®^) for the drop fallspeed in still 

air as in Atlas and Ulbrich (1977) and Q^(D) = CD” for the microwave 
attenuation cross section as in Atlas and Ulbrich (197**). In the latter 
expression C and n would be dependent on radar wavelength and temperature. 

From these definitions it can be shown easily that each of the 
quantities W, R, , A, and can be expressed in terms Z, and m 
through the general form 

where the coetiicient a^, exponent e^^ and function given in 

Table I with F(a) equal to the complete gamma function, ft = a^/r(7+m) 

m m 

and = (AD^) = 3.67+m. The functions FoCD^.m) and F.(D^,m) are defined 

mom no A o 

in the table in terms of the manerical integrals involving v(D) and Q (D) 
in R and A, respectively. Also shown in brackets are the results for ^^gCm) 
and i.i^(m) when the aforementioned power law approximations for v(D) and 
Q^^(D) are used in Fp(D^,m) and F^(D^,m). It is immediately apparent that 
the Ujj for W, F, and are all dependent solely on m so that for these 

quantities the effects of variation in size distribution shape (i.e., 
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variations in m or G) are contained completely in the The same remarks 

apply to the v<|^(m) and li^(m) in the case where the power law approximations 
are used. These can therefore be used to adapt the rain parameter 
diagram to distributions other than exponential (m = 0). 

In making such an adaptation, the ordinate (or Z isopleths) and 

isopleths are considered invariant and the labels on the remaining 

isopleths are multiplied by the ratio U jj(m)/Pjj(0) appropriate to the 

rainfall parameter and size distribution in question. These ratios for W, 

^ and are shown in Fig. 2 over the range -2 ^ m < 2. Also shown is 

the variation of G with m and t.ie ratio V j^(m)/Pjj(0) for the power law 

approximation to v(D). As an example of how to use these ratios, to 

construct a rain parameter diagram for m= 2, the Z and D lines are left 

o 

unchanged and the W, I , N.p, and R labels are multiplied by the ratios 
)jj^(2)/vi^(0) found from Fig. 2, viz., 0.76, 1.22, 0.^19, 12.1, and 1.31, 
respectively. This example (m=2) is shown by Ulbrich and Atlas (1982a) to 
correspond to thunderstorm rain so that for this type of rain the W labels 
are reduced by 2'<X, the E labels are increased by 22t, the labels are 
reauced by about 50> and the labels are increased by more than an order 
of magnitude. 

The ratios p^(m)Ai/^(0) which would be used to relabel the microwave 
attenuation isopleths are almost identical to the Up(m)/yp(0) in Fig. 2 
in the case where the power law approximation to Q^(D) is used. In fact, 
for m > “2, ^ < 3*22 cm and -10 T £ 20^C the M^(m)/P^(0) are almost 
independent of X and T, the maximum difference among these ratios being 
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only about 5t. For m = -2 the maxifflum difference is only 14}. 

The accuracy of the use of the power law approximations to v(D) and 

Q^(D) for relabelling the isopleths on the rain parameter diagram has been 

determined by calculating the and using the numerically integrated 

results for Fj^(D^.m) and F^CD^.m) in Table I. The results so obtained are 

functions of both and m and will be displayed in later work. It is 

sufficient here to note that the differences between the results obtained 

by numerical integration and those found using the power law approximations 

are only a few percent over the range 0.02 ^ £ 0.5 cm for m > -2. For m 

= -2 the differences amount to about 10% over the same range of D 's. 

o 

The remaining measurable on the rain parameter diagram, viz., can 

be adjusted for size distribution variations using the results of Ulbrich 
and Atlas (1982a). The latter work shows that this quantity can be 
approximated closely by the form 


Z 


DR 




(5) 


where is in db and is in cm. Only the coefficient is dependent 
on distribution shape and its variation with m is shown in Fig. 2. Since 
Zpp is a function of only for given m, these results can be used to 
transform the scale on the rain parameter diagram to a Z^p scale. 


4. Simulations of Remote Rainfall Measurement 


Raindrop size spectra have been used by many investigators to simulate 


remote radar measurements so as to discover the relationships between these 
measurables and the rainfall parameter of interest. The earliest and best 
known of these simulations uses raindrop size spectra to calculate the 
radar reflectivity factor Z and the rainfall rate R from which data a 
relationship of the form Z r aR^ is found by standard regression methods. 
The accuracy of prediction of such an empirically determined relation can 
then be tested by finding Z from the drop size spectrian, calculating R from 
Z in the empirical Z-R relation, and comparing the result so obtained with 
R found from the drop spectrum. Such a simulation involves a comparison of 
rainfall rates determined using a simulated remote single-measurement Z-R 
technique with those found from a simulated raingage measurement. Of 
course, it does not include complications due to effects such as vertical 
and horizontal winds, differences in size and location of the measurement 
volumes of the radar and raingage, instrument errors in either the radar or 
raingage, and other effects. In other words, the comparison is between two 
presumably perfectly measured quantities and the accuracy of prediction 
which results from such a simulation is the best that could be expected in 
nature with the use of a single-measurement technique. 

Examples of such approaches to rainfall rate measurement simulation 
are far too nuraercus to review here. It is sufficient to show as 
illustrations of the kinds of results obtained from such analyses those 
given by Ulbrich and Atlas (1977), Atlas and Ulbrich (1977), and by Ulbrich 
(1981). Each of these investigations uses raindrop size spectra collected 
with a momentum disdrometer which were analyzed to determine empirical 
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relations between R and each of the remote measurables Z, Z, and A. For 
the latter measurable several radar wavelengths were considered and the 
results are discussed in Atlas and Ulbrich (1977) and reviewed in the next 
section; in this section only those results for the wavelength Xs 3*22 cm 
are illustrated. The Z-R, E-R and A-R relations found from these empirical 
analyses were then tested for accuracy of prediction in the manner 
described above. 

The results of the ccmiparisons of the rainfall “ates calculated from 

the empirical relations with those found from the size spectra 

(Rgctual' shown in Fig. 3- The upper, middle and lower sets of data 

points correspond to found from the empirical Z-R, --R, and A-R 

relations, respectively. It is clear that the accuracies of prediction of 

these three single-measurable empirical relations are comparable. The 

average absolute deviation (AAD) (i.e., the average deviation without 

regard for sign) of from Rgctual relation, 37f for 

the E-R relation, and 24% for the A-R relation. Similar results are 

obtained for other remote measurables and rainfall parameters, i.e,, the 

accuracy of prediction of empirical relations between any of Z,5' or A and 

any rainfall parameter, such as R, W, D , etc. is similar. Only in the 

o 

special case of X 1 cm does the empirical A-R relation produce 
significantly better agreement between calculated and actual rainfall 
rates. The fortuitous circumstances which produv-e this good agreement are 
described by Atlas and Ulbrich (1977) and reviewed in the next section. 

The origins of the limitations on the accuracv of prediction of these 
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empirical relations have been discussed earlier in Section 2 where it was 

concluded that it requires at least two remote roeasurables to determine 

rainfall parameters accurately. Support for the latter concl iSion has been 

presented by Ulbrich and Atlas (1975, 1977, 1982a) and Ulbritn (1981) who 

use drop size spectra to simulate dual>roeasurement techniques. The pairs 

of measurables which have been considered in these investigations are (Z, 

A), (Z, I) and (Z, Z^^). The first of these pairs consists of those 

measurables determined by a dual-wavelength radar with one attenuating 

wavelength and a second for v^ich the Rayleigh approximation is valid. 

Alternatively, a radiometer may be used to estimate A albeit a 

path-averaged value (see Section 5). The second corresponds to a system 

employing a radar and a lidar, laser or optical transmissometer . The third 

pair represents those measurables acquired by a dual-polarization radar. 

In simulations of these dual -measurement methods, the size distribution is 

assumed to depend on two parameters (such as and in ’'q. (2) with m 

given) so that each of the remote measurables and the rainfall rate are all 

functions of these two parameters. The experimental size spectra are used 

to compute values of the two measurables in question which are then 

combined by simple elimination between the two relevant defining equations 

(see Ulbrich and Atlas, 1978) to find the two distribution parameters. 

These distribution parameters are then used to find a calculated value of 

rainfall rate (R „, ) as deduced from the pair of measurables. The result 
calc 

so obtained can then be compared with the actual rainfall rate 

found from the drop size spectrian in the same fashion as for the empirical 
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analyses described above. 

Comparisons of rainfall rates formed by such dual measurement 
simulations with the actual values are shown in Fig. 4. The upper set of 
data points corresponds to a (Z,A) simulation described by Ulbrich (1981) 
with the attenuating wavelength X = 1.25 cm. The AAD of the from the 

Rgctual only 5H for this set of data. The second set of data points 
also corres'Tonds to a (Z.A) method but with X = 3.22 cm and the AAD for 
this case is also 5%, indicating that the accuracy of this dual-measurement 
method is not dependent on the choice of the attenuating wavelength. 
Comparison of the accuracy of this method with that obtained above with the 
use of Z or A alone in an empirical relation shows that the improvement in 
accuracy is a factor 7 or 5 relative to the Z-R or A-R empirical methods, 
respectively. 

The third set of data points in Fig. represents a comparison of 

R with for a dual-measurement method employing Z and . Here 

the AAD = 8X so that the improvements in accuracy are factors of ^ and 3 

relative to the Z-R and y.-R methods. A closer inspection of this set of 

data points shows that there is a systematic offset of the R , below the 

calc 

^actual* Atlas (1977) show that this offset is due to 

deviations of the experimental size distribution from the assumed 
exponential form ( i .e . , Eq . (<?) with m = 0). If a gamma distribution is 
used with m - c then the average offset of these data is reduced to zero 
and the AAD becomes less than indicating that the improvement in 
accuracy relative to empirical methods is at least as good as the (Z.A) 
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technique when the correct form of the size distribution is used. The 
reason that the dual-measurement (Z,E) method is more sensitive to size 
distribution variations than tne (Z,A) method is that those distribution- 
dependent factors involving Z in the solution for are more strongly 

affected by such variations than those vdiich involve A. This can be shown 
easily from the equations defining R, Z, A, and Z in Table 1. 

The last set of data points in Fig. 4 corresponds to a simulated 
measurement technique involving Z and Z^^ where the former quantity in this 
case refers to the reflectivity factor for horizontal polarization. This 
simulation is the same as that described by Ulbrich and Atlas (1982) in 
another paper in this volume vrtio find the AAD = 27t when it is assumed that 
m = 0 (i.e., an exponential distribution). However, when it is assumed 
that m s 2 the comparison of R^^^^ with Rgctual 

as indicated by Ulbrich and Atlas (1982a), for this case the AAD s 13X, 
indicating that this dual-measurement method is even more sensitive to size 
distribution variations than those discussed above. That this is the case 
is demonstrated by Ulbrich and Atlas (1982a) using an approach similar to 
that described for the (Z,I) measurement technique using the defining 
equations in Table I. 

It may be concluded from these simulations that combinations of remote 
measurables involving Z and A, such as those used in dual-wavelength radar 
techniques, have potential for producing accurate estimates of rainfall 
rate without the requirement of apriori knowledge of the size distribution 
shape. In other words, for these methods only two remote measurables are 
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needed and it is sufficient to assume that the size distribution is 
exponential. For other techniques such as (Z.i:) and (Z, Z^^^) variatitons of 
the distribution from exponential can produce significant effects which 
should be accounted for if these techniques are to realize their fullest 
potential . 

Ulbrich and Atlas (1982a) describe methods by which the shape of the 
distribution can be determined; one of these involves the introduction of a 
third remote measurable. By so doing each of the parameters N^, D^, and m 
can be determined simultaneously for each data point, thereby allowing for 
changes in size distribution shape from moment to moment in the rainfall. 
One such method might use Z, Z^j^, and T. as the triad of remote measurables 
and the results of a simulation for this case are shown in Fig. 5. From 
these data it is obvious that by accounting for changes in distribution 
shape it is possible to achieve very high potential measurement accuracy. 

Of course, a system involving three simultaneously determined remote 
measurables poses profound problems for implementation. A simpler appro ‘^h 
would be to estimate distribution shape in rainshafts using path-averaged 
measurements as described by Ulbrich and Atlas (1982a) and in a later 
section of this work. Alternatively, use may be made of the empirical 
observation that there is a close relationship between the parameters 
and m within a given rain and between rainfall types which may be written 
approximately as N = 0.2 exp (3.5 where = (AD„)„ = 3.67 m. Ibis 
simple relationship, which is obeyed very well both within a given 
rainstorm and between different rainfall types, may be sufficient to allow 
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for changes In distribution shape without tt.e need for the introduction of 
a third remote measurable. 

The results of these simulations of multi -parameter methods indicate 
that these methods have high potential for determining rainfall rates 
accurately. However, we emphasize that these simulations embody the same 
assumptions as those for single-measurement simulations involving empirical 
relations. That is, it is assumed here that each of the measurables in the 
multiple-parameter scheme is measured perfectly, i.e., without instrumental 
error, and that the effects of complications such as vertical winds, 
turbulence, eti. ., are minimal. As in the case of single-measurement 
simulations, the results presented here represent the ultimate accuracy 
attainable by methods employing more than one remote measurable. In spite 
of these simplifying assumptions these simulations reveal the origins of 
the failure of empirical methods and demonstrate the potential improvement 
which is possible. 

5. Remote Measurement Methods 

Ihis section deals with the various remote measurement methods which 
exploit radar, microwave links, and/or radiometry either singly or in 
combination. In the case of radar we assume that the reader is familiar 
with or has access to the literature for the accurate measurer:<ent of echo 
power and the determination of the reflectivity factor for single 
wavelength relations to rainfall (Atlas, 1964; Battan, 1973). Dual 
wavelength radar methods generally involve the measurement of echo power at 
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two or more range gates at two wavelengths aimed at the measurement of 
rainfall attenuation. Dual polarization radar or the so called 
differential reflectivity technique provides yet another approach to obtain 
a parameter which is sensitive to a representative drop size. 

Rainfall can also be deduced from attenuation measured by either a 
radar or a continuous wave transmitter viewing one or more retroflectors 
through the rain, or simply transmitting on a one way path to a receiver. 
Such one way paths may be from ground to ground, to aircraft, or to space. 
Also, microwave radiometers viewing the sun or the absorbing medium itself 
provide excellent proxy measurements of attenuation, and when combined with 
radar or other means of deducing the rain path-length, give average rain 
rates. And of course, we have also alluded to the utility of optical 
extinction in previous sections. 

All of these methods and combinations thereof will be discussed in 
this section along with examples of experiments which demonstrate their 
utility and limitations. Because of their simplicity and the abundance of 
the literature on the subject, we begin Mth a brief discussion of 
radiometric techniques followed by methods of measuring path integrated 
attenuation or deducing it from cither or both microwave radiometers and 
radar . 

5. 1 Radiometry 

For meteorological applications, radiometry was originally introduced 
as a proxy for the path integrated attenuation. Since it c«nes into the 
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present work in a variety of contexts, we shall discuss it briefly here. 

The following treatment follows that of Hogg and Chu (1975). Radiometers 
tracking the sun or simply viewing the rain filled sky were introduced 
largely by the space communications community to obtain proxy data for 
attenuation. We shall not treat the suntracker because it is too 
restricted in its pointing requirements, although it has a considerably 
wider dynamic range than passive radiometers in measuring attenuation and 
provides an excellent calibration for the latter. 

The brightness temperature, T^, measured by a passive microwave 
radiometer may be expressed as 

Tp = Ta [1 - exp (- Ox)] (6) 

where is the apparent temperature of the absorbing rain medium and a its 
attenuation coefficient in cm~^ ; x is distance in cm. Recall that 

olcm"’) = 2.3 x 10"^ A (db km"S (7) 

Thus the total attenu:^tion is 

AL = -10 Log (1 - T^/Tg) (8) 


where L is the total path length in km. 
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Hogg and Chu discuss various methods of obtaining the apparent medium 
temperature T . Probably the best approaches are to use either a 

oi 

suntracking radiometer or a space earth communications link. These provide 
measures of true attenuation (l.e., AL) which, with measured values of 
permit estimates of T , Despite the fact that the temperature of the rain 

a 

varies with height and season, measurements ranging from North Carolina to 

Ottawa gives values of T between 270 and 273 K (Strickland, 197**; 

8 

Ippolito, 1971). This is due to the scattering of the cold sky by the 
drops, an effect which increases with both rain rate and radio frequency 
(Zavody, 1974). Hogg and Chu conclude, that for statistical purposes, one 
may use T = 273K for conversion of T. to attenuation at frequencies of the 

8 D 

order of 20 GHz over much of North America. For individual measurements of 

path attenuation, one needs to use more accurate approximations to T . 

Measurements of attenuation deduced frcmi a 16 GHz radiometer are 

compared to simultaneous 15.3 GHz attenuation values as measured with the 

ATS-5 satellite in Fig. 8 (Penzias, 1970). The accuracy is excellent 

although the maximum is only 5 db. Similar measurements by Strickland 

(1974) discussed below show somewhat greater variance. 

It is clear from Eq. (8) that as the rain attenuation increases T^^ 

approaches T exponentially. Thus when exp (-cx) s 0.1 (and AL = -lOdb) 

8 

T./T = 0.9; i.e., the T./T curve is near saturation. Vfhile measurements 
D a D 8 

of can generally be made to accuracies of about 0.1 to I^K, the possible 
variations in specifying T^ are such that there is general agreement that 
passive microwave measurements should not be used for attenuations larger 
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than 10 db. To obtain a sufficiently large dynamic range in radiometric 
measurements of rainfall will require two or more wavelengths. Each 
doubling of the specific attenuation coefficient attained by increasing the 
wavelength will result in a twofold increase in the maximum RL product 
where R is average rain rate and L is the total path through rain. Fig. 9 
(after Lu (1978) presents curves of versus R for a 3. 2 cm radiometer for 
various path lengths, L. The calculations include the effects of oxygen 
and water vapor. Note that the curves ell saturate at RL product of about 
1000 (ran h”^ km) and that the useful dynamic range is just about two 
decades in rain rate. Of course, the shorter wavelengths saturate at 
smaller RL values and have correspondingly smaller dynamic ranges. In nine 
comparisons to average rain measured by a surface network Lu (1978) found 
that the radiometer underestimated the observed rain by a mean of only 
about 8X. The rain path length was determined by radar. 

It is useful to note that at wavelengths at which the attenuation 
coefficient is close to a linear function of rain rate, we have found (not 
shown here) that the measured T^^ is essentially independent of the 
distribution of R along the path. 

5. 2 Path Integrated Attenuation 
5.2.1 Ground to Ground 

Because it was recognized long ago (Gunn and East, 1954) that the 
attenuation of microwaves of about 1 cm wavelength was essentially linear 
with rain rate, (k)llis (1961) and Godard (1965) attempted to measure 
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rainfall from the 0.86 cm attenuation of echoes between two fixed targets 
spaced some distance apart. The idea is that the measurement is then 
independent of the absolute calibration. Harrold (1967) following the lead 
of Coll is and Godard did the most careful study of this approach using an 
0.86 cm radar, two corner reflectors 6.9 km apart, and 8 rain gages along 
the intervening path. He found a mean attenuation coefficient of 0.27 R 
db/kiti in excellent accord with theory. Ratios of attenuation deduced 
cumulative rainfall to gage values fell within ♦ 20> in all but 6 of 23 
cases. 

The radio science literature is replete with measurements of 
attenuation along one way paths from transmitter to receiver at a broad 
range of wavelengths. In an early study comparing theory and measurement 
Medhurst (1965) reported disappointing results. He found very wide scatter 
and a clear tendency for measured attenuations to exceed maximum 
theoretical values (i.e., those corresponding to monodisperse drop sizes 
having the maximum extinction cross-section). In view of the findings of 
Atlas and Ulbrich (1977) showing excellent correlations between attenuation 
and rain rate computed from drop size spectra, and their conclusion that 
the attenuation near 0.9 cm wavelength was very insensitive to the drop 
size distribution, it is at first surprising that the path measurements 
should show such wide scatter. It is also curious that Harrold (1967) and 
his predecessors failed to get better rainfall measurements using the 0.86 
cm radar target attenuation method. Ttie only 3*^ GHz path-average 
measurements which agree well with theory are those of Norbury and White 
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(1972) which were made over a very abort two way path of only 448 m (224 m 
doubled) and used four fast response rain gages along the path giving an 
average spacing of only 56 m between gages. Perhaps more important in the 
context of our subsequent discussion of sampling errors and the effects of 
up> and downdrafts is the fact that 3 of the 4 gages were on 3 n high 
poles, Just 2 m below the microwave path, and they were sampled at 10 sec 
intervals. It is also pertinent to note that the polarization was 
vertical. Their results given iii Fig. 6, show remarkably small scatter and 
agree reasonably well with theory. They found A = 0.24 R compared to A = 
0.20 R db km~^ for vertical polarization reported by Atlas and Ulbrich 
(1977). On the other hand, the data of Semplak and Turrin (1969) at 18.5 
GHz (Fig. 7) show much greater scatter and are more typical of the results 
reported in Medhurst's (1965) review. 

Hogg and (^u (1975) provide the following reasonable explanations for 
the observed discrepancies and scatter. 

Rainfall is apparently overestimateo from attenuation when; 
o Rain cells or maxima pass between gages, 
o Gages have slow response or are sampled infrequently, 
o Gages do not catch all rain; the error increases with windspeed. 
o An updraft occurs within the microwave path thus reducing the mean 
drop fall speed and enhancing the liquid water content and 
attenuation in the path. 

o Oblate drops attenuate horizontal polarization more than vertical; 
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this was not recognized in early experiments, 
o Failure to account for wetting of antenna radomes. 

Rainfall is underestimated when: 

o Downdrafts in the path increase mean drop fall speed. 

0 Polarization is vertical. 

It is noteworthy that Semplak and Turrin (1969) had clear evidence of 
strong convergence and an updraft on October 25, 1967, when they measured 
consistently larger than theoretical attenuations at 18.5 GHz. The data 
for this case are shown by the crosses in Fig. 7. An updraft of only 1 m 
s”^ would explain this apparent discrepancy. Joss et al . (197^), using a 
vertically pointing Doppler radar, also found updrafts when measured 
attenuation aloft exceeded that calculated from DSD at the ground. This 
also explains the factor of 2 excess attenuation (in decibels) found by 
Anderson et al . (19^7) in a well designed experiment at 24 GHz on a 1.95 km 
path in orographic rain along the windward slope of Hawaii near Hilo. We 
also note that recent studies of so-called storm "downbursts” by Fujita 
(1981) provide persuasive evidence that strong downdrafts can occur in 
storms at only a few tens of meters above the ground. Similarly, with 
intense localized convergence, we will find strong updrafts. 

The effects of up- and downdrafts noted above and of turbulence in 
general raises serious questions about the meaning of surface rain gage 
measurements. The rate at which the rain mass reaches the surface locally 
may be either enhanced by downdrafts or reduced by updrafts. Point to 
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point agreement between either attenuation or radar reflectivity above the 
ground and surface gage rainfall is therefore not to be expected except in 
very steady stratiform rain. (See the discussion of the work of Joss et 
al., 197^, in Section 5). However, over some scale approaching the size of 
the storm cell itself, one should expect these effects to average out. 

This suggests that excessive rainfall along one attenuation path should be 
balanced by a deficit elsewhere. 

Several important implications arise from this discussion: (1) 

Provided one accounts for the non-meteorological errors listed above, 
including sampling, there is no reason to doubt that microwave attenuation 
is indeed an excellent measure of the rain mass, and in the absence of 
vertical air motions, of rain rate. (2) However, if one is to use 
attenuation as a measure of local rain rate, which is only physically 
meaningful near the ground in any case, one must also measure the 
difference in vertical air motion between Ine heights of the attenuation 
path and the gages. (3) One must exert caution in utilizing surface rain 
gage statistics to infer fade statistics especially in convective storms; 
attenuation data derived from radiometers or radars elevated along the path 
of the microwave link should generally be satisfactory with the provisos 
noted below. 

In shof't, we are convinced that once one accounts for the various 
error sources noted above, utilizes both vertical and horizontal 
polarization, and measures vertical air motions in the path, we should find 
the excellent relations between attenuation and path averaged rain rate 
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simulated by Atlas and Ulbrlch (1978), and demonstrated experimentally by 
Norbury and White (1972). Indeed, it would appear that a folded 0.86 cm 
microwave path such as that of Norbury and White or one which traverses a 
square of perhaps 100 m on a side and is elevated only 5 to 10 m above the 
ground would provide a quantum improvement on any rain gage system covering 
a comparable area. 


5.2.2 Ground to Ground; Optical Extinction 

The use of electromagnetic radiation in the optical wavelength band 
for the measurement of rainfall has not been used to any great extent 
because of the obvious imitations of the method due to the presence with 
rain of fog and/or cloud droplets. The latter are usually abundant in 
widespread stratiform rain and severely attenuate a light beam. However, 
in heavy isolated showers their effects are not so pronounced and a limited 
number of experiments have been performed in these conditions using 
coherent optical radiation. 

Atlas (1953) was the first to relate optical extinction to rainfall 

rate. Using the Marshall and Palmer (1948) drop size distribution he 

found a theoretical I-R relation of the form T.= 0.31 whereas his 

0 68 

empirical analysis of raindrop size spectra yields E = 0.20 R * . The 

latter results agree fairly well with the empirical analysis of Ulbrich and 
Atlas (1977) who find E= 0.18 

Wilson and P^nzias (1966) measured attenuBvion of a 10.6 pm maser beam 
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by showers along a 2.6 km path in New Jersey and found maximum extinctions 
of 12 db km”^ in rainfall rates of about 50 mm h”\ This agn es very well 
with the theoretical prediction of Atlas (1953)> 

A carefully performed experiment using laser radiation hao been 
reported by Chu and Hogg (1968) in heavy showers. Their theoretical 
relation for X= 0.63 ym was calculated using the Laws and Parsons (19^3) 
size spectra and can be represented very well by the form I = 0.23 
This result agrees well with Atlas (1953) after correcting for the 
difference between the Marshall and Palmer (19^8) and Laws and Parsons 
(19^13) spectra. Chu and Hogg (1968) also present theoretical results which 
include corrections for beam broadening and forward scatter and which can 
be represented by the form Z = 0. 16 Most of the experimental 

measurements of Chu and Hogg (1968) are in better agreement with the 
empirical relation of Atlas (1953) given above than with the latter 
theoretical result. 

The most successful use of coherent optical radiation to measure 
path-average rainfall rate has been reported by Wang et al. (1979). They 
use a system which detects the Airy diffraction-interference pattern 
produced by a raindrop falling through a collimated laser beam with path 
length of 140 m and diameter of 20 cm. Since the system produces a signal 
for every raindrop which falls through the beam that is directly related to 
the size of the drop, the method is not subject to errors induced by the 
assumption of a specific form of the size distribution. In fact, the path 
average size distribution of the drops can be determined from the data. 
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Wang et al. (1979) display results for N(D) for steady rain which are in 
very good agreement with the Marshall and Palmer (1948) distribution. The 
results found from these distributions for the rainfall rate R are compared 
with those deduced from a tipping bucket raingage immediately below the 
laser beam and are shown to be in very good agreement. The system has 
several advantages over other remote measurement systems, the major one 
being that prior knowledge of the size distribution is not required. 
However, its usefulness is restricted to situations where fog and cloud do 
not seriously affect the measurements. 

5.2.3 Earth to Space 

In recent years there have been a great many studies of path 
integrated attenuation and its relation to rainfall by the space 
communications community. Much of the relevant literature is summarized in 
a fine NASA handbook by Kaul et al . (1980), and in a series of 
comprehensive review papers by Oguchi (1981), Fedi (1981), Brussaard 
(1981), Olsen (1981), and Cox (1981). The previously cited paper by Hogg 
and Chu (1975) is also most valuable. An extensive series of papers 
reporting on propagation through precipitation is also found in a European 
Space Agency report on ATS-6 propagation Experiments in Europe (European 
Space Agency, 1977). From our point of view, the most interesting of the 
above papers are those which compare attenuation to rainfall as measured by 
rain gages under the path, to attenuation and or rainfall as deduced by 
radar, and to attenuation as inferred from microwave radiometry. Space 
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permits discussion of only a few of the key references. 

5.2.M Multi-Parameter Measurements in Space Communications 

The importance of reliable microwave earth-space links at frequencies 
in excess of 10 GHz has led to the need for path statistics on attenuation 
and rainfall. This in turn has led to the widespread use of radiometers 
and radars to deduce path attenuation. Strickland (1974) used 15.3 GHz 
radiometer deduced attenuations for direct comparison to 15.3 GHz 
attenuation on a link to ATS-5. Typical results are shown in Fig. 10. 
Standard deviations of the differences range from 0.25 db at 3 db loss to 
1.2 db at 8 db loss. The 1:1 line corresponds to an assumed effective 
medium temperature of 278K; the bounding curves to a variation of + 10K. 

All the data fall nicely within these bounds. Strickland also used a 3 GHz 
colocated radar directly along the same path and computed the 15.3 GHz 
attenuation through the combination of the radar relationship Z = 200 
and the appropriate A-R relation. After adjustment of the radar constant, 
agreement with measured attenuations was said to be "generally good for 
most meteorological conditions." 

A well designed radar/attenuation experiment was also executed by 
McCormick (1972). He measured path attenuation between three transmitters 
mounted on an aircraft flying along a circle of constant radius from the 
receiver site. A colocated 2.9 GHz radar was used to measure reflectivity 
versus azimuth at the elevation angle corresponding to the attenuation 
path, and to compute the attenuation based upon a method similar to that of 
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Strickland (1974). Fig. 11a shows a combined display of the radar PPI 
contoured in reflectivity according to the inset legend, and an azimuthal 
plot of 15.3 (<Hz attenuation in decibels along the periphery for September 
24, 1968 at times shown along the outer scale. Fig. 11b presents 
comparisons of the measured attenuations at 4, 8, and 15 GHz (heavy lines) 
to those calculated from the radar data (thin lines). McCormick found that 
the radar data could be used to predict attenuation well provided the 
hydrometeors are liquid. However, gross overestimates of attenuation 
occurred when either hail or a bright band occurred. 

Goldhirsh (1976) improved upon the aforementioned studies by using a 
well calibrated high resolution radar and a network of rain gages and 
disdrometers under the earth-space path which comprised 13 and 18 GHz 
transmitters on an uplink to receivers on board ATS-6. When a bright band 
was present on the radar, the attenuation integral based upon the radar 
data was terminated Just below the band. In convective storms without a 
bright band, the Integrals were carried along the earth-space path to zero 
echo. Correlations were generally excellent as shown in Fig. 12. 

Goldhirsh also computed Z-R and A-R relations based upon drop size spectra 
measured at the transmitter site, and also upoi DSD's from Marshall and 
Palmer (1948 land Jess et al . (1968). For 18 GHz, calculated attenuations 
differed from measured values by only 1.3 db rras using his own DSD's, and 
by 1.6 and 2.1 db for M-P and Joss et al . , respectively. 

We conclude from the multi-parameter path integrated studies that the 
radar reflectivities and attenuation, either as measured by earth-space (or 
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earth-aircraft) microwave links, or as deduced from passive radiometers are 
highly consistent with one-another when the scatterers/absorbers are radar 
detectable liquid drops. Indeed, the correlations are remarkably good and 
are free of the tremendous scatter found in either attenuation-rain gage or 
reflectivity-rain gage data. This is clearly due to the fact that in the 
path comparisons (i.e., radar, microwave link, radiometer) the measurements 
ore colocated and not subject to all the errors previously discussed in 
using spaced surface gages to represent the rain mass along a continuous 
path some distance above. 

It is useful to note that McCormick (1972) actually used the 
difference between the radar calculated and the measured attenuation as a 
means of calibrating the radar. Since the difficulties of using rain gages 
and disdrometers to calibrate radar are not unlike the aforementioned 
problems associated with their use along microwave paths, the latter method 
commends itself for purposes of radar calibration. 

5. 3 Some Multi-Parameter Experiments in Radar Meteorology 

Until recently few Institutions working in radar meteorology have had 
either the resources or the motivation to conduct multi-parameter 
experiments corresponding to those in space communications. 

It was largely as a result of the impetus provided by the differential 
reflectivity (ZDR) technique (Seliga and Bringi, 1976) that serious 
interest in the enhanced measurement of rainfall by the combination of 
parameters was revived. It iJ therefore not surprising that actual 


experiments of this kind were sparse until the last few years. 

5. 3. 1 Dual Wavelength Methods 

One of the first discussions of the use of a combination of 
reflectivity factor at one wavelength with relative attenuation at another 
wavelength to deduce raindrop sizes was presented by Austin (19^7). Atlas 
(1954) suggested a similar method which would use dual-wavelength radar to 
measure cloud liquid water content. Several of the assumptions made have 
been carried over into much uf the subsequent work: (a) for a range gate 
containing liquid hydrometeors, the reflectivity factor is nearly 
independent of wavelength, (b) the appropriate measure for a 
dual-wavelength radar is over a path interval rather than a single range 
bin, (c) th® attenuation can be related to the liquid water content (or 
rain rate) by means of an empirical relationship. 

One of the primary interests of this method has beer, its potential for 
hail detection (Eccles, 1975; Carbone, 1972). Although earlier work (Atlas 
and Ludlum, ?961) had considered the logaritKn of the ratio of return 
powers from S-band and X-band radar, Eccles and Atlas (1973) showed that 
the derivative of this quantity with respect to range is a more reliable 
indicator of the leading and trailing edges of a hail shaft. However, 
Srivastava and Jameson n977) showed that there are serious problei.;s of 
interpretation of the hail t'gnai for certain combinations of hail size 
I'istributlon parameters. An important consequence of the technique of 
Eccles and Atlas (19V :>) is that in the absence of hail, the interval 
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attenution can be estimated. More detailed descriptions and error analyses 
of the method have been given by Kostarev and Ciernikov (1968)# 

Sulakvelidze and Dadali (1968), Eccles and Mueller (1971), and Eccles 
(1979). 


For rain rate determination, the essentials of the method can be 
deduced from tiie radar equation and a ratio of powers. The return power P 
from a range at wavelength is approximated by 


P(X;^ ) « • * ~ jj-O.f 4 J Aj: (^) WxJ 

J » 


(9) 


where C. is the calibration constant, Z, in mm /m'^, is the reflectivity 
factor and in db km”\ Is the attenuation coefficient at . For a dual 
wavelength radar wiuh X ^ \ ^ and at two ranges r^, r^^ with rj > r^^ a 

power ratio Fp can be measured, where 


yt __ *j ) 

^ T(\„ r.) 

which, from the radar equation, can be expressed as 
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then 
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The interval wera^ed rain rate R or liquid water content W then can be 
estimated by means of empirical laws relating A to rain rate and liquid 
water content. In most applications the longer wavelength, is chosen 
to be non-attenuating. 

As compared with the standard Z-R method, this type of dual wavelength 
algorithm represents an improvement in that the estimate of R is 
independent of the radar calibration constant. Furthermore, the estimate 
is generally less sensitive to fluctuations in the drop size distribution 
(DSD), a reflection of the relative insensitivity in the A-R law to 
variations in the DSD (Atlas and Ulbrich, 1977). Nevertheless, there are 
limitations of the method. Perhaps the most severe problem is one of 
dynamic range. For an S-band/X-band combination, and at low rain rates, 
the method must be applied over fairly long path intervals to attain 
sufficient sensitivity. For an S-band/K-band pair, the major difficulty is 
saturation at the shorter wavelength whenever long path lengths or intense 
rain rates are encountered. 

A second difficulty is the large numbers of independent samples that 
are required to reduce the variance in the power ratio measurement (Eccles 
and Mueller, 1971; Srivastava and Carbone, 1971). The errors can be 


decreased, however, by interval averaging (Eccles and Mueller, 1971). 
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Eccles (1979) has also described a fairly complex procedure t which uses 
averaging in both range and azimuth, and which chooses certain minimum 
points in Z r)/Z (X^, r) through which a smooth curve is fitted. One 
other possibility is to employ a broad band transmitted pulse (Krehbiel and 
Brook, 1979) to increase the effective number of independent samples. This 
strategy, however, has the disadvantage of increasing the receiver noise 
and further restricting the dynamic range. A less severe error source 
arises from the assumption that the reflectivity factor is independent of 
wavelength: the bias that is introduced, however, will become larger as 

the separation between the wavelengths increases or as the interval rain 
rate becomes more variable. 

Despite these problems, the results reported by Berjulev and Kostarev 
(197^) and Eccles (1979) show that this method can achieve more accurate 
rain rate determination than the standard Z-R method if carefully chosen 
averaging schemes are used. 

In the Eccles and Mueller (1971) and Kostarev and Chernikov (1968) 
approach, R is computed from A through an A-R law. In addition, however, 
there is the possibility of first estimating a kind of mean two parameter 
DSD from their equations. This is outlined in the following development. 
Using the foregoing results and the definition of the attenuation 
coefficient, then for an attenuating wavelength, X^ 3 ^^ non-attenuating 
wavelength, X^, pair 


r 


IJ, = -/ \ M <lx » - « f(p/\ 


(14) 
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Moreover, if the radar equation is integrated between and Tj and the 
definition of Z is used then 

5 - 1 

''x ' *•* 


where 


4 

H - ~ — ^ ,IKlio.i3 is the refractivity factor for water, 


p 

backscattering cross section of a drop of 
diameter D in cm, N (D,x) , in cm“\ is the DSD written as a function of 
range and Qj.(X^,D), in cm^. is the total cross section. The limits on the 
D integration are some assumed values and 

If the DSD is taken to be of the form 


Nfp^x)« AJpW 


( 16 ) 


where A is independent of range over the interval £ fj, then a ratio 

of the above equations is independent of (y) and can be solved for 
Writing I(^) (-AD)dD then in the i'^^ range bin, k < i < j, 

can be found from the equation 
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( 17 ) 


Assuming a drop velocity law, v(D); then the rain rate at the i'^^ bin 
becomes 


=r CTT J (18) 

The main drawbacks of this method are the overly restrictive 
assumption on the form of thx, ixiJ and the dependence of the estimated R on 
errors in the calibration constant. Since the interval attenuation is 
found via the Eccles and Mueller (1971) method, the error sources in A are 
the same as previously noted. 

Over an interval of uniform rain rate, (Joldhirsh and Katz (197^) were 
able to show that the two parameters of a mean DSD could be deduced for 
both attenuating/non-attenuating and dual attenuating wavelength 
combinations. To outline the latter method, it can be shown from the radar 
equation and the uniform rain rate assumption that 

i Is* 

= 1 - — £ (19) 

ft 

where 


= F(\ 




( 20 ) 


Notice that from the same set of measurables used in the Eccles and Mueller 
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(1971) approach a different equation is derived. This is a consequence of 
the uniform rain rate assumption. The paraneter 

A can be found from this equation; is obtained from the equation 


- - 

/V, = ;; 

(iTj 

P 


( 21 ) 


where the subscript i can be set equal to either 1 or 2. 

Error analyses of the method have been performed (Goldhirsh, 1975; 
Stogryn, 1975; Wexler, 1976). Essentially the dynamic range of the method 
is no less restrictive than the Eccles and Mueller approach. Moreover, the 
technique is limited to intervals along which the rain rate is nearly 
uniform. This is a more stringent requirement than the Eccles and Mueller 
assumption that the reflectivity factor is insensitive to wavelength. 
Nevertheless, this dual-wavelength technique offers the potential of 
measuring a two parameter DSD which is independent of the calibration 
constant. In this respect it appears to be unique among the 
dual-measurement techniques that have been proposed. 

One of the first reports of experimental results on the dual 
wavelength method is by Kostarev and Chernikov (1968) who used wavelengths 
of 3 cm and 0.86 cm. The rain rate algorithm is essentially the same as 
that derived by Eccles and Mueller (1971). Comparisons were made between 
the average rain rates over a 2.88 km path as deduced by their dual 
wavelength method and by rain gages spaced 250 m apart. The authors 
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reported good agreement between the two sets of data when the rain rate is 
homogeneous over the interval. Inhomogeneous rain rates, however, often 
led to significant errors. The error was explained partly on the basis of 
unmatched antenna ieamwidths; they reported improvements when the beams 
were matched in azimuth. Of course, other errors as discussed earlier must 
also have been operative. 

Berjulev and Kostarev (1974) have reported additional measurements 
made by a 3,2 cm and 0.86 cm radar system. The path-integrated 
attenuations were compared to the same quantity deduced from the raingage 
data. For the 12 storms measured in 1969 the deviation between the 
measurements vary from 7.71 to 43. 8J with an RMS error of 19. 5X. 

Yamada et al. (1978) have designed a dual-wavelength radar system that 
operates at 5 GHz and 14 GHz ( \^ s 6 cm, = 2.14 cm). Although their 
primary interest was in compiling cumulative statistics of the total path 
attenuation, they have also shown several individual graphs of the 
cumulative path attenuation (CPA) versus the range r. Despite the fact 
that the CPA must be a monotonically increasing function of range, the data 
indicate both small and very marked departures from this behavior. They 
note that the difference is probably not due to the unmatched antenna 
beamwidths and suggest possible effects of the frequency dependent 
reflectivities or rain depolarization. To these we should also add the 
possible influence of errors due to statistical fluctuations and, to 
account for the marked decrease in total attenuation, the possible presence 
of large, non-spherical drops or small wet hail (Eccles, 1979). 
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The results of Eccles (1979) have been referred to earlier. It is 
worth noting, especially in light of the regions of ’’negative" attenuations 
reported by Yamada et al. (1978), that the averaging and curve fitting 
technique of Eccles has been designed for the purpose of detecting iiail 
(thereby identifying spurious regions of "negative" attenuation) and 
eliminating some of the statistical errors caused by insufficient sampling. 
The rather poor sensitivity of the S-band/X-band system employed restricts 
the method to cases of fairly long path intervals and to regions of 
moderate to heavy rain rates. For the twelve most intense storms which 
were analyzed, the rain rates deduced from the dual-wavelength method were 
in much better agreement with the rain gage estimates than those derived 
from a Z-R law using a single S-band radar. Despite these improvements, 
the discrepancies between the dual-wavelength and raingage estimates of 
path-averaged rainfall rate are significant. However, as noted in Section 
5.2 and by Eccles, the gage predictions of rain rate over a path can be 
quite poor, especially for a sparsely distributed network. Tlius, the 
technique is probably better than it appears. 

Undoubtedly the most carefully conducted dual wavelength experiment 
performed to date is that of Joss et al . (197^). They used vertically 
pointing 0.86 and 5.6 cm radars, the latter with Doppler capability, five 
surface disdrometers and three rain gages all within a circle of radius of 
about 50 m. Extreme care was taken to average echoes and obtain best 
estimates of reflectivities from which to compute the 0.86 cm attenuation. 
Both radars were also calibrated by comparing the reflectivities at the 


lowest range gate (580 m) with those computed from the surface DSD for the 
entire 2 1/2 hr teat period. The differences in Z at 0.86 and 5.6 cm due 
to variations in DSD were within ♦ 1 db, but should be correlated in 
height, leaving a net expected error in differential attenuation of ^ 0.5 
db km”\ which is reduced further with time integration. Using the DSD 
data and raindrop radar and attenuation cross-sections, they calculated the 
relation Z = 1360 [Z (ram^ m”^), A (db km~^)] as shown in Fig. 13. 

The data points in Fig. 13 show that the measurements generally 
overestimated the DSD computed attenuation coefficients, the more so the 
smaller the reflectivity. This difference was attributed to contributions 
to attenuation but not to reflectivity by liquid water in cloud form. 

(short wavelength radiometer observations should be similarly affected.) 
They also found Doppler evidence of updrafts of 2.3 ms~\ thus leading to 
enhanced attenuation aloft over that deduced from the surface DSD and gage 
data. This is consistent with our earlier explanations for the excessive 
attenuations measured by Semplak a.id TUrrin (1969) and Anderson et al. 
(1947). 

Important ancillary findings of Joss et al . (1974) relate to the 
sampling problems. They found that correlations between the radar measured 
log Z aloft and that computed from the ground DSD decreased sharply with 
increasing height as shown in Fig. 14. Even when adjusting for the lag at 
which the correlations were maximized, the standard deviation of the 
differences was ‘i db between the surface and 1.22 km height. Without lag 
adjustment, it was 7 db. They also found an rms difference of 13il (i.e.. 
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1 db) between DSD Z's calculated from DSD data collected with disdrometers 
spaced only 100 m apart. This emphasizes what we have indicated earlier: 
that even modest displacement in space and time can lead to great 
differences between measurements idiich are expected to be well correlated. 
Indeed , the errors associated with spaced disdrometers and rain gages are 
usually sufficiently great to attribute much of the spread in scattergrams 
to those alleged "ground truth" instruments. 

While Joss et al. (1974) did not use attenuation and reflectivity to 
deduce rainfall, their results are nevertheless highly relevant to this 
problem. In particular, reflectivity, attenuation, and rainfall are well 
correlated when proper account is taken of the sampling problems. But 
cloud water can affect 0.86 cm attenuation (and radiometry) at low rain 
rates and radar reflectivities. 

The most comprehensive multi-parameter measurements made recently are 
those of Masuko et al. (1981). They used two scanning radars and 
radiometers operatings at wavelengths of 0.86 and 3.2 cm. The entire 
facility was installed in an airplane and flown over the ocean in trials of 
a simulated spaceborne rain measuring system. The antennas were either 
pointed toward the nadir or scanned over an angle of ♦ 23° perpendicular to 
the aircraft path. 

In one mode, they compared the attenuation deduced from the 0.86 cm 
nadir reflection from the ocean by assuming a known surface reflection 
coefficient. The attenuation was then converted to rain rate assuming 
established relations. The 3.2 cm radar was used simultaneously to measure 
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the average Z below the bright band. Fig. 15 presents the data and the 

1 

resulting relationship Z : 126 R • which is quite reasonable despite 
the slightly low coefficient. The downward looking radiometers measured 
brightness temperatures as shown in Fig. 16, also compared to the 
attenuation-deduced rain rates as described above. Of course, no rain 
gages were available on the ocean surface, but it is clear that the 
relations found are also reasonable. The correlation coefficient between 
the rain rate deduced from the 0.86 cm attenuation and the 3*2 cm 
brightness temperature is 0.91. However, the correlation with the 0.86 cm 
brightness temperature is only 0.67, probably because the shorter 
wavelength radiometer is affected by non-radar detectable cloud liquid 
water. Of course, at vertical incidence in stratiform rain it is not 
unreasonable that the total path attenuation from the 0°C level down should 
be consistent with the radar reflectivity below the 0®C level since the 
rain is probably vertically homogeneous. On the other hand, errors might 
be expected at the lower rain rates because of the excess attenuation 
through the bright band. In any case, the results show promise for the 
three methods (i.e., radar, path attenuation, and radiometry) although the 
latter two are proxies for one another. 

5. 3. 2 Dual Polarization Methods 

In 1976 Seliga and Bringi proposed a dual polarization technique from 
which a two parameter exponential drop size distribution could be obtained. 
The suggested measuring instrument was an incoherent S-band radar with the 
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capability of measuring the co-polarized return powers along the 
horizontal, P^, and vertical, P^, directions. The method has since been 
used both in estimating the DSD and in discriminating rain from hail (e.g.. 
Cherry et al., 1980); in this discussion only the former will be dealt 
with . 

It should be noted that it is probably Barge (1972, 1974) who should 
be credited with the concepts which triggered the work of Seliga and &"ingi 
(1976) for it was he who first showed that the cancellation depolarization 
ratio (CDR) was correlated with rain rate. He correctly attributed this to 
the increasing number of large oblate raindrops accompanying increasing 
rain rates. However, Humphries (1974) found that the CDR was contaminated 
by depolarization in the intervening rainfall. It was then that Seliga and 
Bringi conceived of using orthogonal linear polarizations from which 
sufficient information could be obtained to deduce the DSD parameters. 

The basic idea of the technique is that the ratio is independent 

of the multiplicative factor of the DSD; therefore the slope A of the 

exponent can be determined. N is then deduced from either P. or P . 

o n V 

Explicitly, from the radar equation and from the definition of Z, 

■p (r.) = J1£L (22) 

% r- J 1 

2 

where '^^(D) in cm is the backscattering cross section of a drop whose 
volume equals that of a sphere of diameter D. The subscript q = {h,v} 
refers to an incident and backscattered field polarized along the q 
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direction. To write in the form above requires the use of several 
assumptions that are discussed below. 

Using Eq. (1) with A = 3.67/D^ it follows that Zjjj^ = 10 log 

where 


E± ^ 


M e*^(-A D)Jp 
j <r;^fp)«x^(-Ab) V2) 


(23) 


P 

From knowledge of C^/C^ and ‘^q(D), one can solve for A and since the 
right hand side of Eq. (23) i;4 a function of D only. The behavior of Z~,- 
and Z^/N^ as functions of is shown in Fig. 17. After determining D^, 
is then computed from the equation 


A/ 

N = (24) 

HC, f aj(») tr^f(-AP) </J> 

In contrast to the dual wavelength method which can be applied only 
over an interval (usually longer than the range resolution), potentially 
the dual polarization method can be used at each range bin. An obvious 
advantage of a method vdiich provides an estimate of the DSD from measured 
quantities is that it does not require an empirical law to relate radar and 
meteorological quantities. For the determination of liquid water content 
or rain rate the success of the method is largely dependent on how 
accurately the DSD can be determined. 


Seliga and Bring! (1978) listed a number of possible error sources in 


49 


the methods Since then many experimental studies and error analyses have 

b.en carried out. What follows is a brief discussion of the major errors 

in the method and a summary of some of the relevant work. 

It is evident from the manner in which A and N are determined that 

o 

both and either C|^ or must be assumed. Since a single antenna is 

used to measure and P^, generally can be determined more 

accurately than or alone. t')f course, the errors in these quantities 

can be minimized by standard calibration methods, by in-situ calibration 

using rain gages and disdrometers (Browning, 1978) and by procedures for 

calibrating the : io log (Z^/Z^) measurements (Seiiga et al., 1981). 

Another possibility is to find another measurable that is independent of 

the aosolute calibration constant. Generally, such measurement schemes are 

not (judsi polarization methods in the normal sense of the word. 

Since the method offers two measurements, only a two parameter DSD can 

be specified. Moreover, the and D^^^^ limits of integration must be 

assumed. Cherry et al (1979), Hall et al. ( i980a) and Uibrich and Atlas 

(1982a) have studied the effect of on the measurement, showing 

that the latter is sensitive to wiien the median equivalent voliane 

max ^ 

diameter is greater than about 0.2 cm. dowever, if D /D > 2.5 then the 

max o 

dependence of Z on D is essentially independent of variations in D 
iJlbrich and Atlas have also shown that if the three parameters of a gamma 
DSD could be obtained, the improvement in the method would be substantial. 
To obtain the additional measurement required, the use of path integrated 
microwave or optical attenuation has been suggested (see Section >i) . 
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It was noted earlier that several assumptions must be made in order to 

express the reflectivity factor as the product of integrated 

over the equivalent volume diameter D. For a more general assemblage of 

particles, this single integral would be replaced by integrations over 

density, size, shape and orientation angles of the particles. To reduce 

this multiple integration to a single integration it is first assumed that 

the rain drops are oblate spheroids with a one to one correspondence 

between D and the ratio of the minor to major axis of the particles 

(Pruppacher and Beerd, 1970; Green, 1975). Secondly, it is assumed that 

the fluctuations of the angle of orientation of the minor axis about the 

vertical are sufflcieiitly small that the rain drops may be assumed to be 

aligned along a common direction (Seliga and Bring!, 1978). This preferred 

orientation of rain drops is, in fact, the main basis of distinguishing 

liquid from solid hydrometeors (Cherry et al., 1980; Hall et al., 1980b). 

Several studies of particle shapes have been made (Bring! and Seliga, 

1977a, 1977b; Seliga and Bring!, 1977). The effects of drop canting angle 

have also been investigated, but to our knowledge none have been used to 

translate such effects into errors in the method. A final error source 

is due to the finite number of independent samples at each polarization. 

To reduce sanspling errors, and P^ should be staggered with a brief time 

interval between them and a much longer time interval between each 

P^, P^ pair. In this way successive P^, P^ measurements are highly 

correlated k^ile each P. , P pair is nearly decorrelated with all 

h V 

others (Bring! et al . , 1978). Under this situation standard errors 

in Z of about 0.075 db can be made with about 50 independent P. or 
UK n 
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Py samples (Hall et al., 1980a). There are also different ways of forming 
the Zjjjj ratio. Bring! et al . (1980) has shown that an average cf the form 
<^DR> ' (where the angular brackets denote a sample mean 

over the received pulses) is asymptotically unbiased and has the smallest 
standard deviation among the estimators which are possible. 

Despite the number of potential error sources in the technique, the 
experimental results have been very encouraging. Basically, four kinds of 
comparisons have been made between meteorological parameters derived from 
the method and the same parameters measured by more direct means. 

Using a relationship for attenuation A, in terms of N and D , the Z»d 

O O UK 

method was used to compute values of N ,D and then values of A at each 

o o 

range bin (Hall et al., 1980c). Summing the estimated values of A along 
the path of an 11.6 GHz satellite downlink provides an estimate of the path 
integrated attenuation with which the directly measured attenuation can be 
compared. For the rather small values of attenuation which were measured, 
agreement between the two sets of attenuation versus time curves is quite 
good, exhibiting a standard deviation of about 0.3 db as shown on the left 
of Fig. 18. In contrast, the attenuation found from an A-Z relation, 
derived by assuming a Marshall-Palmer DSD, was in error by about a factor 
of two. The agreement between directly measured attenuation and that 
deduced from raingage data was even poorer (Fig. 18, right side). 

In an experiment using the CHILL radar (Seliga et al., 1981), rain 
rates derived from the Z^^^^ technique were compared against rain gage 
measurements. The agreement between these sets of rain rates were 
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significantly better than comparisons with the rain rates deduced by means 
of a Z-R law. This held true even vdien a rain gage calibration method was 
used along with the Z-R method. 

Another experiment using the CHILL radar and a disdrometer has been 
reported (Sel’ga et al., 1980). For this case, the temporal record of the 
disd^'oi; 3ter v.?s used to reconstruct the range-profiled rain rate by means 
of a.i estimate cf the storm speed. Again, the Z^p -derived R compared more 
favorably witu the- disdroroeter-derived R than did the R estimated from Z 
alone. 

A more direct comparison between quantities measured by a radar and a 
disdrometer has been given by Goddard et al . (1982). In that experiment, 
the radar measured values of Zj^p = 10 log (Z^/Z^) and Zj^ were compared to 
the values derived from a disdrometer. The conversion of disdrometer data 
to the radar quantities, Z^p, enabled them to dispense with any apriori 
assumption concerning the DSD. They found that the radar measured Z_„ was 

Dn 

on the average 0.3 db smaller than that deduced from the disdrometer, while 
the radar measured Z^ exceeded the disdrometer estimate by 1.6 db. The 
authors note that although the Z^ discrepancy can be explained on the basis 
of disdrometer errors, this is not the case with the Z^p estimates. To 
reduce the bias in radar measured Z^p, the relationship between D and the 
ratio of minor to major axes of the oblate spheroid is modified so that, in 
effect, the smaller drops are taken to be more spherical in shape than 
would be predicted from the equilibrium results of Pruppacher and Beard 
(1970). This modification reduces the Z^p bias to 0.15 db while having a 
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negligible effect on Z^. A partial justification for this procedure is 
found in the experimental data of Jameson and Beard (1982), although the 
data suggest shapes that are more nearly spherical than those of drops in 
equilibrium at larger values of D as well. Other possible errors sources 
are the effects of canting angle, unequal sampling volumes of the radar and 
disdrometer, and the time lag caused by the spatial separation (120 m) 
between the sampling volumes. 

In another recent paper, Stickel and Sellga (1981) report on estimates 
of and W (liquid water content) made by both in-si tu measurements using 
Knollenberg (1976) probes and by means of the method. Excellent 
agreement was found between the two values, but the Zjjg method provided 
a value of W which was about a factor of 10 less than the in-situ 
measurement. Ihe authors explain the source of discrepancy as due to an 
error in the radar calibration constant. 

5. 4 Range-Profiling Algorithms 

There are several cases of interest in which the following situation 
arises: we have available the attenuated radar return powers at each range 

gate as well as an estimate of the total path attenuation (taken at the 
same wavelength and along the same path) and we wish to determine the rain 
rate along the range direction. These are the data that are obtained, for 
example, in the fixed target method and in the radar/radiometer and radar/ 
satellite downlink combinations (see Section 5.2). For the latter two 
sensors it is assumed that the wavelengths are identical. 
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According to the hypothesis, the total path attenuation, AL , can be 
estimated, where 


~ A('x) <lx ( 25 ) 


The path length is equal to vh, where vis the total number of range 
gates and h is the range resolution. 

A second estimate of this quantity can be obtained by means of a 
method first derived by Hitschfeld and Bordan (1954). They have shown that 
for an attenuating wavelength radar, the rain rate at the j bin, R , can 
be expressed as a function of the return powers P (k = 1,.., j), 

, the radar constant for wavelength and the parameters in the A = ^ 

and R = aZ*^ laws. Assuming that a, b, and g are independent of range, 
then we can write this solution in the form 






-1 




• ,'».J 9 •»••• 

1j » J - L 


(26) 


where 




( 27 ) 


where 


^kj ~ I /c#j •as' /or /r-j 

J 
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From the radar equation and the expression for R^, we can infer the 
path- in teg rated attenuation out to the j range bin (J £ v). Ihis is 
simply 

ALj = cok^re Lj zrjA. <28) 

By letting J increase to a second estimate of the total path-integrated 
attenuation is obtained, i.e., 

(AL^y ^ (29) 

where f^ is given by Eq. (27) with j replaced by v. In general (AL^)', 
which is determined from the return powers from the individual range gates, 
is not equal to AL of equation (25) which is found by means of a 
path-integrated measurement. However, by adjusting , or in the 
expression for f^, we can insure the agreement between the two. The 
parameters so determined can then be inserted into Eq. (26) which yields 
the range profiled rain rate. Details of this algorithm are given in 
Meneghini et al. (1982). 

A simple error analysis of this method has been carried out (Meneghini 
et al . , 1982). It was shown that an adjustment of the calibration constant 
leads to range profiled rain rates that are independent of errors in the 
calibration constant and in the Z-R and A-Z laws, but dependent on the A-R 
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law and errors in the path-integrated attenuation measurement. In the 
analysis, however, it was assumed that the R-Z and A-Z relationships were 
constant along the entire path. Although a and ^ were treated as randomly 
varying from path to path, the parameters b and B were assumed to be 
constant and unbiased. 

A somewhat different approach to this problem has been given by Lu 
(1980). Although the method was developed for a radar/radiometer sensor 
pair, it can be applied unchanged to the other measurement systems cited 
above. 

r 8 

From the radar equation and an A =s Z law, we can write 



(30) 

where the exponent of the equation is found from the previous (i-1) 

equations. These equations are now subjected to the constraint that the 

radar path integrated attenuation AL = h (A -»•,.+ A ) be equal to the 

> v-1 

path integrated attenuation derived from the radiometric intensity 
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measurement, t . That is, 

P 

=i.(A,* ■■■■■* A^„) ( 31 ) 

This constraint allows us to vary either ^ or g in the above equations 
until this equality is satisfied. This being accomplished, the sets of 
{A^} and {Z^} quantities follow immediately. 

It should be noted that the first to propose a technique of this type 
were Hitschfeld and Bordan (1954). They suggested that the parameters in 
the rain rate estimate could be adjusted so that agreement was obtained 
between this and a rain gage reading deep in the storm. This method serves 
to bound the errors in the rain rates determined up to that point. In 
essence, the only difference between their method and the others here 
described is that a different measurable is used as a constraint. 

Finally, it should be mentioned that methods of this type may be of 
some use in the context of the dual attenuating/non-attenuating wavelength 
method. Since the interval attenuation is monotonically increasing 
function of the interval length, it follows that the fractional standard 
deviation of this quantity will generally decrease with increasing path 
length. For a variety of situations, this interval may be comprised of 
many adjacent range gates. Under these circunstances , it is not difficult 
to show that with minor changes in the range profiling algorithms described 
above, the rain rates can be estimated at each range gate within the 


interval of interest. 
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6. Some Inferential Methods 

While it would be desirable to develop a remote measurement scheme 
i^ich could be completely automated, we recognize that at least for the 
time being, human intervention can often greatly enhance the accuracy of 
measurements. For example, the use of a Range Height Indicator (RHI) 
display provides a wealth of information concerning the horizontal and 
vertical homogeneity of the rainfall, the vertical gradients of 
reflectivity due to growth, evaporation, or wind shear drop size-sorting, 
or as a result of displacement of the rain streamers out of the plane of 
observation due to cross winds. All of these features relate to both 
errors in sampling at the surface and to ways of correcting surface 
measurements for such effects. 

It is also clear that in the case of stratiform rainfall, considerable 
information about the nature of the DSD and precipitation rate is contained 
in the reflectivity profile across the bright band (BB). Until now, this 
information has not been exploited because we have lacked adequate vertical 
resolution to view the BB from the side. However, theoretical studies of 
the evolution of the DSD and the accompanying Doppler spectra across the 
melting layer by Ekpenyong and Srivastava (1970) suggest that the depth of 
the BB is sensitive to both the median volune particle size as well as to 
the precipitation rate, assuming no complexities associated with 
aggregation or breakup. While the absolute depth of the BB also depends 
upon the temperature lapse rate within the melting layer, the depth of the 
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BB may be defined by the layer encompassed within the two -6 db points 
above and below the height of the maximum reflectivity. While it is 
generally acknowledged that the depth of the BB increases with 
precipitation rate, at least to rates of about 10 nin/hr (which is very high 
for stratiform precipitation in any case), we know of no attempts to use 
that as a measure because of the inadequate resolution of most radars. 
However, it should not be overlooked at least as a qualitative indicator 
with high resolution radars or at vertical incidence where high resolution 
may be achieved with short pulse lengths. 

Another BB indicator of precipitation intensity which must be 
considered qualitative for the time being is the ratio of the peak 
reflectivity of the BB, ^p, to that above, and/or below, The ratio 

Pp/pg should increase with precipitation rate because the rate of snow 
crystal aggregation increases with the precipitation rate. Similarly, the 
greater the aggregation, the greater the breakup, and thus the larger the 
ratio Hp/n p. While these relationships have been noted qualitatively, we 
are unaware of any quantitative studies which demonstrate them clearly. 

Many years ago it was also noted (Newell et al . , 1957) that the 
circular depolarization ratio (CDR) was often much smaller in the melting 
layer than in either the snow above or the rain below. This was attributed 
to both the flatness of the snowflake aggregates and their tendency toward 
horizontal orientation as anticipated by Atlas et al . (1953). With the 
advent of the dual polarization differential reflectivity (ZDR) method, 
similar results have again been reported in the BB (Hall et al . , 1980b). 
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It is also pertinent to note that Cherry et al. (1979) have noted that 
the ZDR peaks at a height a few hundred meters below the peak of the BB; 
indeed, it appears to occur near the height where the reflectivity itself 
‘’as attained the value corresponding to rain. They attribute this to the 
large differential polarization effects of oriented wet flakes much as 
described above. But, in order that the peak of the BB differ so much from 
that of ZDR, it seems clear that the particle size weighting of the 
reflectivity must differ markedly from that for ZDR. The suggestion is, 
therefore, that the peak ZDR occurs essentially just above the level at 
which the largest particles in the size spectrum melt completely. If this 
is correct, then the difference in height between the peaks in Z and ZDR 
should be a measure of the breadth of the particle size distribution. 
Evidence in support of this contention should be attainable from high 
resolution vertically pointing Dtjppler measurements. 

It would thus appear that the combination of the reflectivity ratios 
across the bright band and ZDR measurements should also be at least a 
qualitative indicator of aggregation, breakup, size spectrum breadth, and 
of rain intensity. This can readily be determined by correlating either 
CDR and ZDR in the BB to Z and or rain rate below; similar correlations 
should be sought with the depth of the layer between the peak Z and ZDR. 

Another feature which should be explored is the effect of size sorting 
due to wind shear on the drop size spectrum, Z, and ZDR. Gunn and Marshall 
( 1955 ) have modeled and Atlas and Plank (1953) have observed the history of 
the DSD, reflectivity, and rain rate due to wind shear. When precipitation 
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generating cells are sufficiently far apart, it is well known that the 
initial drops to reach the surface at the leading edge of the cell are the 
largest and thus produce Z*s larger than expected with unsorted DSD's of 
the same rain rate; thr converse is true at the trailing edge. Effects 
such as these account for much of the scatter about the mean Z»R regression 
relationship, and indeed, also for systematic shifts in the Z-R relation on 
the rain a paraneter diagram. Examination of high resolution RHI displays 
of Z and ZDR by the first of us (private communication with S. M. Cherry, 
1981) has shown evidence of drop size sorting. Because much of the theory 
on which the ZDR method is based assinnes an exponential DSD, observations 
such as these should permit the determination of when such an assumption is 
valid and the development of correction algorithms for use when it is not. 

Finally, while this work has concentrated entirely on rainfall, a few 
words are in order concerning the discrimination of rain from hail. It is 
apparent from the RHI profiles of reflectivity factor and differential 
reflectivity presented by Hall et al. (1980b) that the combination of these 
two measurables can be used to distinguish rain from hail or water-coated 
ice particles. Although a quantitative, automatic method for identifying 
hydrometeoi type using these measurables has not yet been established, the 
method shows promise. This is particularly true when the measurements are 
combined with subjective criteria by the observer. 

Another multi-parameter method for identification of hail is the dual- 
wavelength hail detection scheme of Eccles and Atlas (1973) referred to in 
Section 5. The usefulness of the technique has been investigated by 
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Srlvastava and Jameson (1977) who suggest that the ratio of the 
reflectivity factors at two wavelengths is a better measure o ' hydrometeor 
type. The latter method has been further investigated by Ullrich and Atlas 
(1982b) and implemented for Colorado hailstorms by Jameson and Heymsfield 
(1980). This method also not yet been developed fully so that automatic 
identification of hydrometeor type and size can be made. 

7. Summary and Conclusions 

In this paper we have attempted to review the rationale behind the use 
of multi-parameter observations aimed at rainfall measurements of improved 
accuracy. We simulate the potential improvements attainable with various 
combinations of remote measurables, discuss the various concepts proposed, 
and review those experiments which have been attempted for these and 
related purposes. 

In Section 2 we summarize t. e long history of Z-R and drop size 
distribution (DSD) measurements which finally led to the conclusions that 
radar rainfall relationships suffered from both systematic variations and 
scatter due to deviations in the drop size and number from well established 
DSD's such as that due to Marshall and Palmer (19**8). To overcome this 
problem Atlas and Qimela (1957) developed a rain parmneter diagram which 
required the specification of two parameters. Ulbrich and Atlas (1978) 
extended that diagram by adding several other parameters including 
microwave and optical attenuation, but requiring the DSD to be exponential. 
The use of the latter rain diagram by Ulbrich and Atlas (1982a) showed that 
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the accuracy of rain rate specification could be improved significantly by 
the use of a third variable which accounted for the breadth of the 
distribution. This was done through the use of a gamma DSO. Accordingly, 
in the present paper (Section 3), we generalize the Ulbrich-Atlas rain 
parameter diagram still further. 

In Section 4 we use DSD's to compute or simulate the remotely 

measurable parameters reflectivity, Z, specific attenuation. A, optical 

extinction, ^ , and differential reflectivity, Z and the rain rate, R, 

UK 

and median volume diameter, D , We then demonstrate that R calculated from 

o 

pairs of remotely measured parameters such as (Z, A) or (Z, ^ ) decreases 
the errors by factors of about 6 and 4, respectively. Finally, the use of 
three parameters (Z, produces virtually perfect agreement. While 

these simulations indicate the considerable potential value of the use of 
multiple parameters, they fail to simulate the various sources of error in 
the actual measurement systems. 

The remainder of the paper is devoted to a review of the var. ' US 
remote measurement concepts and illustrative experiments which demonstrate 
their capabilities and limitations. We fopus first on path-average 
microvove attenuation. Despite the disappointing agreement between most 
experimental attenuation measurements ard the theoretical values expected 
from surface rain gage data which disturbed Medhurst (1965), we agree with 
Hogg and Chu (1975) that most of the errors are due to rain gage sampling 
problems and to vertical air motions which produce differences between the 
rain rates in the elevated paths and at the level of the gages. When these 
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problems are eliminated as in the well designed 35 GHz experiment of 
Norbury and White (1972), excellent results are found, so much so, that wc 
believe that a similar system modified to cover an area of some 100 m on a 
side would be far superior to any presently conceivable raingage system. 

We also emphasise that the sampling and vertical air motiovi problems 
which have plagued the above experiments also account in large part fur 
much of the discrepancies in radar-rainfall experiments. In short, there 
is considerable doubt about the "truth” of rain gages. 

It is also shown that passive microwave radiometry is a fine proxy for 
path integrated attenuation up to values of about 10 db. It has the great 
advantage of not requiring either a reflector or receiver at the end of the 
path, thus allowing measurements from a single site, or even with the same 
antenna as may be used for simultaneous radar measurements. However, in 
order to achieve sufficient dynamic range one must use a multiplicity of 
wavelengths. Radiometric measurements of path average rainfall by Lu 
(1980) show good agreement with gage values. 

The availability of stable microwave beacons on communications 
satellites has opened the door to their use in path attenuation 
measurements. Moreover, the need for rainfall and attenuation statistics 
for earth-space communications has given new impetus to the development of 
a variety of schemes to obtain proxy attenuation data from weather radars 
and radiometers. Various Investigators have used radar-rainfall and 
attenuation-rainfall relations to transform radar measurements to path 
attenuatxon. Very good agreement has been found in this way except when 
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the radar beam traverses the bright band or hail. Agreement is improved 
further when the basic relationships are adjusted for the actual DSD below 
the path. Similarly good agreement has been found between radiometer- 
deduced and actual attenuation. 

In the area of optical extinction surprisingly little has been done 
since the early predictions of Atlas (1953). One very well done experiment 
by Chu and Hogg (1968) finds good agreement between their measurements and 
the empirical findings of Atlas (1953). The disadvantages are that it is 
restricted to short paths and, in the case of light rain, may be affected 
by fog, dust, and/or pollution. 

We then proceed to the variety of dual wavelength and dual 

polarization (2^^) techniques which have been treated in radar meteorology. 

For a variety of reasons including excessive signal fluctuations, 

difference .'.ri reflectivity and/or unmatched beams at the two wavelengths, 

non-Rayleigh scatterers such as hail, and the modest attenuation at some 

wavelengths and small rain rates, the methods have not produced very 

promising results over short paths. However, Eccies (1979) appears to have 

achieved reasonable success in using the method over a larger range azimuth 

2 

domain in comparison to cumulative rainfall measured by a 1600 km surface 
network of gages. In essence, however, we believe that the method reduces 
to a path averaged measurements technique. A variant on the above is the 
method of Goldhirsh and Katz (197^) which assumes a uniform rainfall rate 
over some path to arrive at the effective exponential DSD. The requirement 
for a uniform rain rate is surely too stringent except in steady stratiform 
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precipitation. 

Wiile the work of Joss et al. (197^) was not aimed at evaluating the 
dual wavelength method, it did demonstrate that reflectivity, attenuation, 
and rainfall rate are highly correlated when proper account is taken of the 
sampling problems. But undetected cloud water can produce excessive 

i 

attenuation at short wavelengths at small rain rate and reflectivity. Ihe | 

most exciting recent multi-parameter experiment is that of Masuko et al. 

(1981) who used radars and radiometers at 0.86 and 3.2 cm in an airborne 1 

laboratory. Using the 0.86 cm radar in a nadir reflection mode from the | 

sea, they measured attenuation. The rainfall deduced thereby was well 
correlated with the 3.2 cm reflectivity below the melting level and with 
the radiometric brightness temperatures at both wavelengths. 

IXial wavelength techniques and radiometric methods are more promising 

‘i 

When they are used to estimate attenuation. The latter is then used as a * 

i 

constraint on the retrieval of range profiled rain rates from the radar 
measurements. The best results in gate-by-gate measurements have been 
achieved with dual polarization or differential reflectivity (Z^^). 

However ’ven these have failed to meet their full potential because 
rainfall often does not behave according to the apriori assumptions. An 
accompanying paper (Ulbrich and Atlas, 1982a) shows that the use of a third 
remote parameter in addition to Z and Z^^^^ offers great promise. 

In many cases, accuracy can be greatly enhanced and ambiguities 
resolved by personal examination of the radar displays which depict the 
nature of the bright band, cell spacing, homogeneity, and drop size sorting 
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One very important field of multi-parameter measurement techniques 
which has not been covered in depth in this review is that which includes 
polarization diversity methods other than those discussed in this work in 
connection with The use of polarization methods for the measurement 

of precipitation is covered in a paper by Humphries and Barge (1982) to be 
presented at this symposium. The methods are not sensitive to canting 
methods and have therefore been used primarily as a means of identifying 
hydrometeor type. However, recent work which has been directed toward 
relating polarization diversity effects to rainfall parameters and other 
precipitation characteristics indicates that the methods have considerable 
promise. 

With a growing appreciation of the needs and the capabilities of the 
various techniques and growing evidence of the willingness to invest the 
necessary effort and resources to do the job right, the outlook for highly 
improved remote rainfall measurements seems bright. 
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TABLE OF SYMBOLS 


Units 


Definition 


(db km" ) 

(db km'^ 

(db km 

(db km'M 

. 6-b ,b “3. 

(mm h m ) 

(see Table 1) 


r 2-n. 
(cm ) 


5 -6 

(w m mm ) 

r s -6. 
(w in mm ) 


( S -6. 
(u in mm ) 


!^db cm 


microwave attenuation. 

microwave attenuation at wavelength 

X. and distance x. 

1 

microwave attenuation at wavelength A. 

path-average microwave attenuation. 

coefficient in Z-R relation ^ . 

coefficient in relation between para- 
meter X, Z and D (Table 1) . 
o 

average absolute deviation, 
total path attenuation, 
bright band. 

exponent in 2-R relation Z = aR^. 

coefficient ir power law approximation 
to total attenuation cross section 
Q^(D) = con. 

radar calibration constant. 

radar calibration constant for hori- 
zontal polarization. 

radar calibration constant for vertical 
polari zation . 

coefficient in power law approximation 

- = r n 

m o 

circular depolarization ratio. 

equivalent spherical raindrop diameter. 

size distribution parameter median volume 
diameter . 

upper limit of raindrop diameters in drop 
size distribution. 

lower limit of raindrop diameters in drop 
s i ze distribution . 

drop size distribution. 

exponent in relation between parameter 

X. : and D (Table 1). 
o 

numerical integral in definition of A 
(Table 1). 
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Symbol Units 


Definition 




numerical integral in definition of R 
(Table 1). 


f. 

3 


dimensionless function relating rainfall 
rate to power return to radar. 


H 

(cm^) 

h 

(km) 

1(A) 


iKl^ 


L 

(km) 

L 

V 

(km) 

m 


N{D) 

(m cm 

N(I),x) 

(m cm 

N 

o 

(m ^ cm 

N 

o 

(m cm 

>< 

o 

(m cm 

N . 

Ol 

(m cm 

Nf 


n 


r(A..rJ 

(watts) 


(watts) 

P (r.) 

(watts] 


( watts ) 


dimensionless third moment of the rain- 
drop mass distribution 


d {] D*'N(D)dD 


N(D) dP 


H - 



2 


range resolution; interval between range 
gates . 


1(A) = 


o(D) cxp(-AD) dD. 


D 


= 0.93, refractivity factor for water. 


total path length. 

path length spanned by v range gates. 

exponent in gamma size distribution func- 
tion N(D) = exp(-AD) . 


raindrop size distribution. 

raindrop size distribution at distance x 
from radar. 

raindrop size distribution parameter. 


path- or interval -average size distribution 
parameter. 

rcindrop size distribution parameter at 
distance x from radar. 

raindrop size distribution paramete?* in 
i th range gate from radar. 

total raindrop number concentration. 

exponent in power law approxim«Tr ion to ^ 
total attenuation cross section Q^(D) = CD . 

power returned to rr- ! at wavelength X. 
from range r. . 

power return to radar from range r. at 
horizontal polari ation. ' 

power return to i adar from range r. at 
vortical polari . at ion . 

average power return to radar at hori- 
zontal polari z at ion . 
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Symbol 

Units 

Definition 

<p > 

V 

(watts) 

average power return to radar at vertical 
polarization. 

Q^CD) 

(cm^) 

total attenuation cross section of rain- 
drop of diameter D. 

Q^(D.X) 

(cm^) 

•jtal attenuation cross sec .ion of rain- 
drop of diameter D at wavelength X. 

R 

(mm h”^) 

rainfall rate. 

R. 

1 

(mm h~^) 

rainfall rate at ith range gate from radar. 

R 

(mm h 

path average rainfall rate. 

\ctual 

(mm h'^) 

rainfall rate found from raindrop size 
spectrum data. 

*^calc 

(mm h ^) 

rainfall rate found from empirical rela- 
tion or from multiple-measurement simulation 

r . 
J 

(m) 

radar range. 

T 

a 

ri'' 

J 

physical absorber temperature. 

\ 

(K) 

microwave brightness temperature. 

T. 

1 

2 

(watts km ) 

Ti = P (X., r^)r^^ 

v(D) 

(m s 

raindrop fallspeed in still air. 

W 

(g 

liquid water content. 

W 

(g m'^) 

path- average liquid water content. 

X 

(see Table 1) 

any of the parameters W, R, I, A, N„, N 
(Table 1). 

X 

(cm, km) 

distance. 

z 

(mm*’ m*^) 

radar reflectivity factor. 

Z(X.,r.) 

, 6 -3, 

(mm m ) 

radar reflectivity factor at wavelength 
X^ and at range r from radar. 

^DR 

(db) 

differential radar reflectivity factor. 

<Z > 
DR 

01 

m 

6 

6 

m 

T’ 

1 

P 

(db) 

pulse-average differential reflectivity 
factor. 

(AD ) = 3.67 + m. 

o m 

g 

exponent in A-Z relation A = C Z . 

6 = J/T(7 + m). 

m m 

power ratio 
P(Xj,r.) 

P(Xj,r^) P(X^.T .1 • 
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Symbol 


Units 


Definition 


r(a) 


■Y 


E . . 
IJ 


x.x, 


, -1 -0.67 

(in s cm 


(cm“ ra ) 
(cm“ m 
(cm“ m 
(cm'^) 

(cm'S 

(cm) 


complete gamma function 

OD 

f a-1 -u , 
r(a) = u e du . 

o 

coefficient in power law approximation 

to drop fa 11 speed v(D) = y 

IT nT -0.67 

Y = 17.67 m s cm 

'ij ' ‘ ‘ I* j 

G . . = 0 . 5 i = i . 

11 

radar reflectivity above bright band, 
radar reflectivity below bright band, 
peak radar reflectivity in bright band. 

raindrop size distribution parameter. 

path- or interval -average size distri- 
bution parameter. 

wavelength . 


u^(D^,m) 


V 


r 

(db mm 

z 

(km” S 



0 

(cm ) 


(cm“) 




(cm") 

0 .(D) 

(cm") 


dimensionless pseudomoment of the drop 
size distribution for parameter X. 
(Table 11. 

total number of range gates. 

coefficient in A-Z relation A - t, Z . 

optical extinction . 

absorber attenuation coefficient. 

radar backscattcring cross section of a 
raindrop of diameter D. 

radar backscattering cross section at 
horizontal polarization for a raindrop 
with equivalent diameter D. 

radar backscattering cross section at 
vortical polarization for a .aindrop 
with equivalent diameter D. 


g] 


rill 


Parameter X, coefficient a^, exp< 

X = a p 

X X 

in terms of reflectivity factor Z(mm 


(D ,m) Z D 
o o 

^ m '^) and med: 



B 

m 


a 


m 

r(7 + m) 


F/Oo.m) 


-1 „-(ni . n . 1) 

o 







84 


r.uGii^iNu Frtut Oi_mi-.A i»vi * i- = 


FIGURES 


I 

i. 


A 



85 


ORIGINAL FAGS IS 

OF POOR QUALITY 


Fig. la: Rain paraneter diagram: base diagram. Radar reflectivity 

factor versus rainfall rate with isopleths of median volume 
diameter Dq, size distribution parameter Nq, liquid water 
content W, and optical extinction ^ shown as heavy solid lines, 
light solid lines, heavy dashed curves, and light dashed curves, 
respectively. Also shown as an inset is the differential 
reflectivity factor Z-_ as a function of D„. [from Ulbrich 
and Atlas (1978)3. 
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Fig. 1b: Rain parameter diagran: attenuation overlap 

Isopleths of microwave attenuation A(db km 
T = -10, 0, 10 and 20°C are shown as dashed, 
and heavy solid curves, respectively, [from 
(1978)]. 
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Fig. 1c; Five typi •il empirical Z-R relations plotted on a simplified 

version oi the base diagram in Fig. la. Also shown as a shaded 
area is the part of the diagram within vAiich fall all of the 69 
empirical Z-R relations listed by Battan (1973) [after Ulbrich 
and Atlas (1978)). 
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— 1 5 

Fig. 2: Coefficient C (db cm ) in power law approximation, 

dimensionless third moment of the mass distribution G, and 
ratios of the pseudomemants b^(m)/b (0) for X = H, W. J’ , 

Np and A versus of the exponent m in the gamma drop size 
distribution. Each of the curves is labeled with the function 
to which it corresponds. 
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Fig. 3^ Coro pai iso ns of calculated rainfall rate with actual 

rainfall rate R for single-measurement simulations 

using a set of arop size spectra described by Ulbrich and Atlas 
(1977). The upper, center and lower sets of data points 
correspond to calculated from empirical Z-R, E-R 

relations, respectively, in each of which the remote measurable 
(Z, ^ or A) was found from the size spectra data. In all cases 
Ractuai determined directly from the size spectra. 
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Fig. 4: 


Comparisons of calculated rainfall rate R with actual 
rainfall rate For dual-measurement simulations using 

the same set or drop size spectra as in Fig. 3. The upper, 
second , third and lower sets of data points correspond to R . ^ 
calculated from simulated measurable pairs [Z, A( ^ = 1.25 cm?, 
[Z, A( A = 3.22 cm)], [Z, Z] and [Z, respectively, for 

vAiich all remote measurables were found from the size spectra 


data. In all cases R_q^„ n is determined directly from the 
size spectra. ^ 
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Comparison of calculated rainfall rate actual 

rainfall rate R i for a triple-measurement simulation 
using the same set or drop size spectra as in Fig. 3. 
data points correspond to R calculated from a simulated 
measurable triplet (Z, Z__,f ; in Miich all measurables were 
found from the size spectra. 
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Fig. 6: 


Attenuation versus rainfall rate for ten rain showers during 
the summer of 1971 in Slough, England. The crosses are measured 
data points, the solid curve represents a theoretical 
calculation of attenuation as a function of rainfall rate using 
the Laws and Parsons (19**3) drop size spectra, and the dashed 
curve is a parabola fitted to the data by least squares [from 
Norbury and White (1972)]. 
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Plot of observed total path attenuation y (db) of 

microwaves with wavelength equal to 1.62 cm^vlrsus path-average 
rainfall rate as measured by raingages for 23 rain storms 

in New Jersey during the summer of 1967. The circles and 
crosses are measured data points, the latter corresponding 
to one storm for which there was evidence of an updraft. The 
line composed of one long dash and two short dashes is the 
theoretical calculation of Gunn and East (1954), the long dash- 
short dash line and solid line are power law and linear least 
squares fits to the data, respectively, and the upper and lower 
dashed lines are the maximum and minimun attenuation-rainfall 
rate relations for monodisperse drop size distributions, [from 
Semplak and Turrin (1969)]. 


^MtAS (total path) in decibels 
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Fig. b: Total attenuation measured by a 16 GHz passive radiometer versus 

attenuation of the 15.3 GHz ATS-5 satellite beacon for a rain 
shower in New Jersey. The data points are measured values and 
the dashed line has a slope of 1.1 to account for the difference 
in frequency of the two signals, [after Penzias (1970)]. 
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Fig. 9: Theoretical curves of microwave brightness temperature T^^dc) 

for a wavelength of cm versus path-average rainfall rate 
R(nin h~ ) and total path length L(km). Elevation angle of the 
radiometer is assumed to be 4^. [after Lu (1978)] 
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Fig. 10: Comparison of attenuation at 15.3 GHz calculated from 

radiometric measurements of sky noise temperature with that 
measured directly using the ATS-5 satellite beacon. The solid 
circles represent measured values for a storm in Ottawa, Canada. 
The straight line of unit slope assumes an effective absorber 
temperature of 276K and solid curves above and below that line 
correspond to absorber temperatures of 268K and 288K, 
respectively, [after Strickland (19711)]. 
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Fig. 11: (a) Contours of reflectivity factor Z for a thunderstom 

observed by a 2.9 GHz radar scanning in azinuth at an elevation 
angle of 'v 6.5^. Also shown along the periphery is the total 
attenuation measured using an aircraft beacon at 15 GHz circum- 
navigating the storm. 
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Fig. 11: (b) Comparison of attenuations measured directly (heavy curves) 

at 4, 8. and 15 GHz using aircraft beacons as in part (a) with 
that computed from the radar data (light curves), [after 
McCormick (1972)]. 
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Comparisons of attenuations measured by ATS-6 satellite 
receivers at 13 and 18 GHz (heavy curves) with that predicted 
by radar data combined with raingage and disdrometer data for 
a storm in Virginia during 1975. [from Goldhirsh (1976)]. 
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) 

Fig. 1i: Attenuation A versus reflectivity factor for a storm in 

Sudbury, Massachusetts, during 1971. The open circles, solid 
circles and crosses represent experimental points measured by 
a dual-wavelength radar with wavelengths of 0.86 and 5.6 cm at 
three different altitudes in the storm. The straight line is 
an empirical fit to A-Z data computed from disdrometer data 
at the ground. The curve at the top illustrates the relative 
number of radar profiles used in computing the experimental 
data points [from Joss et al . (1974)]. 
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Fig. 14: Cross correlation functions between the logarithm of 

reflectivity factor measured aloft at several altitudes by 
radar and the logarithm of reflectivity factor calculated 
from disdrometer data at the surface versus time lag t. The 
radar and disdrometer are the same as those in Fig. 13. 

[from Joss et al . (1974)]. 
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Fig. 15: Average reflectivity factor Z measured below the bright band 

by an X-band radar versus rainfall rate deduced from 
attenuation measured by a K^-band radar. Both nadir-pointing 
radars were colocated on an^aircraft flying over the ocean, 
[from Masuko et al. (1981)]. 
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Fig. 16: Microwave brightness temperature at X-band (left figure) and 

at K -band (right figure) versus rainfall rate deduced from 
attenuation of the K -band radar in Fig. 15. The solid circles 
are measured data pofnts. The solid curves represent least 
squares fits of second degree polynomials to the data points, 
[from Masuko et al . (1981)]. 
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Fig. 17 


Theoretical curves of differential reflectivity factor 
ZogCdbl and reflectivity factor at horizontal polarization 
Zy ) divided by NQ(m”'^cm~' ) versus median volume 

diameter for an exponential drop size distribution 
with maximum drop diameter D = 0.8 cm. [from Seliga 
et al. (1981)] 
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Fig. 18: Left: Comparison of total path attenuation deduced from dual- 

polarization data with A_, that measured directly from an 
11.6 GHz beacon on the ESA geostationary OTS. The solid circles 
and crosses are measured data points, the latter corresponding 
to cases where the melting layer contribute at least half the 
total attenuation. The measurements were made in two storms 
during October 1979 in Southern England. 

Right: Comparison of raingage-derived total path attenuation, X'' 
with that measured directly X^. All data were averaged over 
time intervals of 4 minutes. Other details the same as 
at left, [from Hall et al . (1980c)]. 
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Time sequence of (a) Z^p and (b) as determined from radar 
and disdrometer data. Solid curves calculated from 30s 
averages of disdrometer data using backscattering cross 
sections at horizontal and vertical polarizations as described 
by Seliga and Bringi (1979). Solid circles are radar-measured 
30s averages with the antenna fixed. Open circles are radar- 
measured data with the antenna scanning at 1° s~ . The dashed 
curve is the recalculate*.’ using the modified backscattering 
cross sections of Goddard “t" al . (1982). [from Goddard et al . 
(1982)]. 



